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SECTION  I 

INTRODUCTION  AND  SUMMARY 


This  final  report  describes  a study  of  mathematical  algorithms  for  inferring 
atmospheric  temperature  and  humidity  structure  from  remotely-sensed  radiances.  The 
primary  objectives  were  to  determine  the  efficacy  of  these  algorithms  when  applied  to  the 
sensing  channels  of  the  DMSP  E-  and  H-packages,  to  select  those  which  retrieve  maximal 
information  consistent  with  limitations  set  by  instrument  noise,  and  to  implement  the 
selected  algorithms  as  computer  codes  suitable  for  inverting  real  DMSP  data. 

The  study  succeeded  in  defining  the  capabilities  of  the  DMSP  E-  and  H-pack- 
ages, which  proved  to  be  quite  good  with  respect  to  temperature  retrieval,  but  marginal,  at 
best,  with  respect  to  humidity  retrieval. 


Three  of  the  six  humidity  channels  were  found  to  be  redundant,  and  the  remain- 
ing three  too  strongly  absorbing  to  infer  low-level  moisture  except  for  relatively  dry 
atmospheres.  Because  of  this  outcome  it  was  decided  (with  approval  from  the  Contract 
Monitor)  to  de-emphasize  some  of  the  planned  study  tasks  and  adopt  the  additional  task  of 
determining  how  the  humidity  sensing  performance  could  be  improved  by  the  addition  of 
two  window  channels. 

The  added  window  channels  (at  795  cm-1  and  900  cm-’)  were  analyzed  to 
determine  their  effectiveness  in  inferring  the  total  water  vapor  column  thickness  of  the 
atmosphere  and  the  effective  surface  temperature.  When  the  window  channels  are  used  in 
this  way,  the  remaining  three  (non-redundant  DMSP)  humidity  channels  are  required  only  to 
distribute  the  total  water  in  altitude.  It  was  found  that  the  selected  window  channels  are 
too  opaque  to  accurately  infer  total  water  in  most  cases;  consequently,  the  augmented  set 
of  five  channels  yielded  no  more  information  than  the  three  non-redundant  DMSP  channels. 
However,  it  is  likely  that  improved  results  can  be  obtained  by  incorporating  these  same 
channels  as  additional  sounding  channels,  that  is,  by  performing  a simultaneous  five-channel 
inversion,  since  the  window  channel  weighting  functions  are  large  near  the  surface. 
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In  summary,  the  study  provides  a comprehensive  analysis  of  the  remote  temper- 
ature sensing  capabilities  of  the  DMSP  E-  and  H-packages,  and  a computer  code  that  will 
retrieve  the  maximum  information  available  on  vertical  temperature  structure.  The 
principal  results  of  the  study  relative  to  humidity  retrieval  are  an  identification  of  the 
shortcomings  of  the  H-package,  and  suggested  approaches  to  be  followed  in  future  research 
to  improve  the  capabilities  of  DMSP.  Some  of  the  by-products  of  the  study,  for  example 
the  empirical  transmittance  functions  developed  for  the  DMSP  humidity  channels,  are 
useful  in  themselves.  Also,  the  comparisons  of  different  temperature  retrieval  techniques, 
although  based  on  inversions  of  synthetic  DMSP  band  radiances,  are  applicable  in  a quali- 
tative sense  to  all  remote  vertical  temperature  sounders. 

1.1  STUDY  PLAN 

A plan  of  study  was  adopted  early  in  the  contract,  which  defined  the  specific 
tasks  to  be  performed,  the  basis  of  selection  of  candidate  retrieval  algorithms,  and  the  key 
questions  to  be  addressed  in  evaluations  of  the  candidate  algorithms.  The  plan,  in  some 
respects,  was  more  ambitious  than  the  contract  work  statement,  and  was  not  followed  in 
every  detail.  However,  it  is  worth  presenting  the  plan  because  it  defines  the  scope  of  the 
study,  the  context  of  discussions  in  the  following  sections,  and  the  motivation  for  some  of 
analyses  and  comparisons. 

Table  I lists  the  planned  major  tasks  of  the  study.  The  purpose  of  the  first  two 
tasks  is  to  identify  a 2-by-2  matrix  of  specific  candidate  algorithms  judged  to  be  the  best 
representatives  of  linear  and  nonlinear  methods,  used  with  and  without  a priori  statistical 
data.  The  two  algorithms  in  the  "statistical"  class  use  empirical  covariances  to  regularize 
(stabilize)  the  r<**rieval  solutions;  statistical  regression  methods  were  excluded  from 
consideration.  Tasks  3 and  4 compare  the  capabilities  of  the  four  candidates,  as 
demonstrated  in  inversions  of  synthetic  radiances,  and  select  the  "best"  one  as  the  recom- 
mended procedure. 

The  questions  addressed  in  evaluating  and  comparing  the  candidate  retrieval 
techniques  are  listed  in  Table  2.  All  of  the  questions  applicable  to  temperature  inversions 
were  answered,  except  the  one  dealing  with  instrument  noise.  That  is,  we  did  not  add 
pseudo-random  noise  values  to  the  synthetic  radiances  used  in  the  retrievals,  although  we 
chose  values  for  regularization  parameters  which  should  effectively  damp  any  instabilities 
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resulting  from  the  noise  equivalent  radiances  of  the  DMSP  channels.  In  studying  the 
retrieval  capabilities  of  DMSP-H  and  the  two  window  channels,  we  addressed  parts  of  the 
first  four  questions,  but  the  answers  obtained  are  mostly  qualitative. 


Tcfcle  I MAJOR  TASKS  OF  DMSP  RETRIEVAL  ANALYSIS  PROGRAM 


1.  Comprehensive  Testing  of  Known  Retrieval  Techniques 
Applied  to  SSH  Temperature  and  Humidity  Channels 

2.  Selection  of  "Best"  Technique  in  Each  Class  of  Retrieval 
Method 

3.  Simulation  of  Retrievals  via  "Best"  Techniques,  to  Obtain 
Statistical  Measures  of  Accuracy  and  Information  Content 

4.  Definition  of  Recommended  Procedure  for  SSH  Retrieval  of 
Temperature  and  Humidity  Profiles 

5.  Coding  of  Recommended  Retrieval  Algorithms 


Table  2 QUESTIONS  ADDRESSED  IN  THE  STUDY  OF 
SPECIFIC  RETRIEVAL  ALGORITHMS 

1.  Overall  Accuracy 

2.  Maximum  Achievable  Resolution  in  Vertical  Structure  (Tem- 
perature and  Humidity) 

3.  Dependence  of  Solutions  of  Initial  Guesses 

Magnitude,  Variance  and  Shape  of  Residuals 

4.  Effects  of  Statistical  Constraints  (Covariances) 

Effects  of  Using  "Bad"  Statistics 

5.  Effects  of  Noise 

6.  Optimal  Method  for  Selecting  Guess  Humidity  Profile 

Effects  of  Using  Statistics 
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SECTION  2 

STUDY  OF  TEMPERATURE  RETRIEVAL  TECHNIQUES 
2. 1 SELECTION  OF  RETRIEVAL  TECHNIQUES 

We  initially  selected  ten  candidate  techniques  to  be  studied  with  respect  to 
temperature  retrieval.  These  are  listed  in  Table  3.  Note  that  the  techniques  are  grouped 
into  four  classes,  according  to  whether  they  are  linear  or  nonlinear,  and  whether  or  not 
they  use  statistical  data  in  the  regularization  of  solutions.  Each  class  is  represented  by  one 
technique  except  the  nonlinear,  nonstatistical  class.  The  study  first  examined  the 
nonlinear,  nonstatistical  methods  except  the  one  due  to  King,  since  it  was  being  studied 
independently  by  Dr.  King  at  AFGL.  After  extensive  retrievals  of  synthetic  radiances,  we 
selected  the  "best"  of  these  methods,  which  reduced  the  number  of  candidates  to  four;  i.e., 
one  in  each  class.  The  procedure  was  then  repeated  to  determine  the  "best"  of  the  four. 

J.  I.  F.  King's  nonlinear  technique'*  obtains  an  inverse  solution  of  the  radiative 
transfer  equation.  With  this  exception  the  selected  nonlinear  techniques  are  direct 
methods:  they  solve  the  radiative  transfer  equation  in  the  forward  direction,  initially  for  a 
guess  profile,  then  iteratively  for  revised  guesses  determined  by  differences  between  the 
computed  and  measured  radiances.  The  direct  methods  differ  from  each  other  only  in  the 
relaxation  formulas  used  to  revise  the  temperature  profile  at  each  stage  of  iteration.  We 
regard  the  Chow  method  as  a generalized  and  improved  version  of  the  well-known  Chahine 
method,  which  is  why  the  latter  is  not  listed  as  one  of  the  selected  methods. 

The  linear  methods  use  inverse  solutions.  These  are  generally  regularized  least- 
squares  solutions  of  an  equation  set  consisting  of  a linearized  form  of  the  radiative  transfer 
equation  written  out  for  each  of  the  sounding  channels.  The  names  in  the  last  column  of 
Table  3 identify  two  well-known  matrix  solutions:  the  "minimum-information"  solution 
obtained  by  Foster**  and  independently  by  Twomey^,  and  the  more  general  statistical  form 
derived  in  different  fashions  by  Foster  , by  Strand  and  Westwater  , and  by  Rodgers  . 
These  formal  solutions,  however,  can  be  implemented  in  different  ways.  Usually  one  solves 
for  a highly  over-specified  temperature  profile  (defined  at  many  more  levels  than  these  are 
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Table  3 RETRIEVAL  TECHNIQUES  TO  BE  COMPARED 


NON- 

STATISTICAL 


STATISTICAL 


NONLINEAR 

Direct 

SMITH1 

CHOW2 

FLEMING3  (ANALOGS) 

NEW  (MODIFIED  TWOMEY4) 

AND  VARIANTS  (2) 

Inverse 

KING5 


Direct 

FLEMING  STATISTICAL3-* 


LINEAR,  INVERSE 


Regularized  LeabE-5quares 


► FOSTER6.  TWOMEY7 


Statistically-Regularized 

Least-Squares 

FOSTER6,  STRAND  & WESTWATER8 


spectral  bands)  or  for  the  coefficients  in  an  expansion  of  the  temperature  (Planck  function) 
profile  in  empirical  orthogonal  functions;  both  techniques  avoid  preselection  of  pressure 
levels,  which  would  restrict  the  form  of  the  retrieved  profiles.  One  also  can  choose  among 
several  methods  of  linearizing  the  radiative  transfer  equation,  and  whether  or  not  to 
represent  the  transfer  equation  in  perturbation  form  (use  the  method  with  iteration). 
Details  of  how  the  two  linear  retrieval  methods  were  implemented  are  given  below. 


Note  that  the  two  nonlinear  methods  attributed  to  H.  E.  Fler,  ing  in  Table  3 are 
identified  as  analogs  of  the  two  selected  linear  methods.  The  Fleming  methods  use  linear 
relaxation  equations  derived  from  the  requirement  that  the  mean-square  error  in  Planck 
function  change  (as  prescribed  by  the  relaxation  formula)  be  a minimum.  As  a result,  the 
relaxation  formulae  are  simple  algebraic  analogs  of  the  corresponding  matrix  solutions  of 
the  linear  methods;  this  mathematical  relationship  is  noted  on  page  14. 


2.2  DEFINITION  OF  ALGORITHMS 

2.2. 1 Direct  Nonlinear  Methods 


The  direct  nonlinear  methods  start  with  a guess  temperature  profile  which  is 
us^d  to  compute  corresponding  Planck  radiance  profiles  B..,  where  i denotes  the  spectral 
channel  and  j the  atmospheric  pressure  level.  The  radiative  transfer  equation  in  quadrature 
form, 


Bij  W>f 


= 1,2,. 


M, 


is  used  to  obtain  the  radiances  I.  corresponding  to  the  guess  profile.  In  the  present  study,  M 
is  equal  to  6,  but  J can  be  as  large  as  desired.  In  our  retrieval  simulations,  the  atmospheric 
temperature  is  defined  at  100  levels;  the  value  j = J = 101  corresponds  to  the  boundary  term 
of  the  radiative  transfer  equation  (the  surface  temperature). 


In  the  selected  direct  nonlinear  methods,  either  the  difference 


Al.  = I.  (meas)  - lj 
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between  the  measured  and  computed  radiances,  or  their  ratio 


R.  = I.  (meas)/l. 

i i i 

is  substituted  into  a relaxation  equation  to  obtain  an  adjustment  AB..  to  the  computed 
Planck  radiance  profiles  B...  This  gives  M new  Planck  radiance  profiles,  one  for  each  value 
of  i. 


The  next  step  in  the  procedure  is  to  apply  the  inverse  of  Planck's  equation  to 
the  revised  radiance  profiles  to  obtain  M corresponding  temperature  profiles.  Finally,  the 
M different  temperature  profiles  are  averaged  in  some  fashion  to  obtain  a single  tem- 
perature profile.  This  profile  replaces  the  initial  guess,  and  the  entire  procedure  is 
repeated  until  an  appropriate  convergence  criterion  is  satisfied. 

The  following  is  a summary  of  the  different  relaxation  equations  used  by  the 
selected  methods,  and  their  origins.  The  simplest  formula,  due  to  Smith*, 


AB..  = 
•J 


(I) 


increments  the  Planck  radiances  equally  at  all  altitudes.  Twomey^  proposed  a general 
relaxation  formula  for  directly  solving  the  linear  Fredholm  equation;  if  we  modify  it  to  deal 
with  the  radiative  transfer  equation  in  its  actual  nonlinear  form,  the  result  is 


W..  Al. 

•J  i 


max 


:.  rw..i 

J L i]J 


B.. 

'J 


(2) 


where  the  first  factor  in  the  right  hand  side  is  the  usual  weighting  function  W..  normalized 

3 'J 

to  its  maximum  value.  Fleming's  "new  nonlinear"  technique  uses  the  formula 


AB..  = 


W.. 

•J 


W.j2  . a 


Al.; 


a = constant 


(3) 


H 


which  is  the  result  of  requiring  a linear  relaxation  equation  that  minimizes  the  mean-square 
error  in  AB..,  and  making  simplifying  assumptions  about  the  statistics  of  the  atmospheric 
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temperature  profile  (the  same  ones  leading  to  Foster's  minimum  information  solution).  The 
Chow  method  uses  the  formula 

AB..  = (Rj^  - I)  B..;  K = constant  (4) 

For  K = I , this  reduces  to  the  simplest  ratio-type  (Chahine-type)  relaxation  formula 


B..  + AB.. 
•J  iJ 


(4a) 


However,  in  the  Chahine  method  each  j (level)  is  associated  with  a particular  i (channel), 
and  the  number  of  levels  equals  the  number  of  channels. 

The  Smith,  Fleming  and  Chow  nonlinear  methods  all  prescribe  Smith's  temper- 
ature profile  averaging  scheme: 


M /M 


where  T..  are  the  temperature  profiles  obtained  by  solving  the  inverse  Planck  equation  for 
the  new  radiance  profile  estimates  B...  It  seemed  reasonable,  therefore,  to  use  the  same 
scheme  in  the  method  adapted  from  Twomey  . Note  that  approximating  the  radiative 
transfer  equation  by  a linear  Fredholm  equation  is  symbolically  equivalent  to  replacing  B.. 
by  a function  b.  which  depends  only  on  j;  i.e.,  the  case  treated  by  Twomey  does  not  require 
averaging  over  the  spectral  channels. 

Note  that  both  the  Fleming  method  and  our  modified  Twomey  method  (Equa- 
tions 3 and  2),  prescribe  a AB.j  which  has  the  same  vertical  distribution  as  the  weighting 
function  W...  This  would  seem  to  be  a desirable  feature,  since  according  to  the  radiative 
transfer  equation  the  desired  AB„  information  contained  in  a change  Al  is  distributed  as 
Wjj.  That  is,  we  should  vary  B more  at  altitudes  where  a unit  change  will  produce  a large 
change  in  computed  radiance,  than  at  altitudes  where  a unit  change  would  have  a relatively 
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small  effect  on  I.  However,  if  this  weighting  is  used  in  the  relaxation  equation,  should  it  be 
used  again  in  the  averaging  formula?  It  is  conceivable  that  this  "double-weighting"  will 
produce  undesirable  biases  in  the  retrieved  temperature  profile.  Hence,  we  chose  to  also 
study  the  effects  of  using  the  alternate  (unweighted)  temperature  average 

M 

WE  Tu  <« 

i=l 

with  the  Fleming  and  modified  Twomey  methods. 

Table  4 lists  the  combinations  of  relaxation  and  averaging  formulas  used  in 
studying  the  direct,  nonlinear,  nonstatistical  techniques. 

Table  4 ALGORITHMS  FOR  THE  DIRECT  NONLINEAR,  NONSTATISTICAL  METHODS 


NAME 

GRAPHICS 

RELAXATION 

AVERAGING 

SYMBOL 

FORMULA 

FORMULA 

SMITH 

SMH 

Eq.  1 

Eq.  5 

TWOMEY-LIKE 

TWM 

Eq.  2 

Eq.  5 

ALTERNATE  TWOMEY-LIKE 

TW2 

Eq.  2 

Eq.  6 

FLEMING  NEW  NONLINEAR 

FLM 

Eq.  3 

Eq.  5 

NEW 

NEW 

Eq.  3 

Eq.  6 

CHOW 

CHW 

Eq.  4 

Eq.  5 

It  is  interesting  to  note  that  Fleming's  result  (Eq.  3)  without  the  parameter  a 
can  be  derived  in  a very  simple  manners  Apply  the  identity 


AB..  W.. 
|J  iJ 
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to  the  radiative  transfer  equation  in  quadrature-perturbation  form,  and  require  that  the 
factor  AB../W..  equal  a constant.  The  transfer  equation  is  then  identical  to  Equation  3 with 
a = 0.  In  fact,  it  is  questionable  whether  a non-zero  value  of  a (which  has  exactly  the  same 
meaning  as  the  single  regularization  parameter  in  Foster's  minimum  information  solution)  is 
required,  since  it  is  generally  not  necessary  to  regularize  the  direct  nonlinear  methods. 
They  are  inherently  more  stable  because  they  are  direct  and  also  because  the  temperature 
averaging  operation  tends  to  reduce  the  effects  of  noise. 

2.2. 1 . 1 Statistical  Direct  Nonlinear  Method  (Fleming  Statistical 


We  include  a derivation  of  the  relaxation  formula  for  this  case,  since  we  do  not 
know  if  one  exists  in  the  literature.  The  method  seeks  a linear  formula 


AB..  = C..  Al.;  i = I,  2, ...  M;  j = I,  2, . . . J (7) 


giving  the  least-squares  best  fit  to  an  ensemble  of  data  AB..^  = B..^  -<B„>  and  Al.^  = 1.^  - 
< I >.  These  represent  a set  of  actual  (in  situ,  measured)  Planck  function  profiles  and  the 
corresponding  computed  or  measured  upwelling  radiances;  k denotes  the  ensemble  members 
and  < > a mean  over  the  ensemble.  Hence,  it  is  required  to  find  the  matrix  C which 
minimizes  the  error 


i > 

KC).  £ UB.jk  - C..  4|.k); 


Differentiating  S with  respect  to  the  matrix  element  C„  and  setting  the  result  equal  to 
zero  gives 


rv 

iC  2A,ik  (ABijk  " Cij  = 0 
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or 


assuming  the  noise  has  zero  mean.  From  Equations  7,  10,  12,  and  13,  the  relaxation  for- 
mula is 


AB„  = A|. 
|J  ' 


Fleming  proposed  to  use  this  equation  together  with  Smith's  channel  averaging  scheme. 

Note  that  the  Fleming  statistical  relaxation  formula  revises  B..  at  a particular 

frequency  i and  level  j according  to  the  weights  W for  all  levels,  but  the  weight  for  level  n 

is  scaled  in  proportion  to  the  expected  correlation  between  AB.  and  AB...  This  will  tend  to 

in  ij 

introduce  expected  structural  features  to  the  retrieved  profiles,  which  tends  to  increase 
their  average  accuracy  while  sacrificing  the  ability  to  recover  unusual  or  anomalous 
profiles;  of  course,  this  statement  needs  to  be  quantified. 

If  we  have  no  knowledge  of  temperature  correlations  or  the  vertical  distribution 

of  Planck  function  variance  (or  choose  to  reject  the  use  of  a priori  information)  then  the 

2 

matrix  S..^  becomes  the  identity  matrix  times  , the  variance  of  B.j.  Then  Equation  14 
reduces  to  Equation  3 with  1 


y s..  w. 

LmJ  ijn  in 
n=  I 


J J 

EW..  T S"  W-  + i 

ij  ijn  in  E. 

>1  n=l  1 


(14) 


a = 


(15) 


which  may  be  called  the  ratio  of  noise  to  "signal"  power.  This  is  the  minimum  information 
case,  analogous  to  Foster's  minimum  information  linear  inverse  solution. 
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2.2.2 


Linear,  Inverse  Methods 


The  linear  methods  require  that  we  first  put  the  radiative  transfer  equation 


into  a linear  form;  here  Bv  denotes  Planck's  function,  T(p)  the  temperature  profile,  Tv(p) 
the  transmittance  of  the  atmosphere  above  level  p,  and  x is  any  convenient  single-valued 
function  of  pressure.  We  will  use  the  linearization  scheme  suggested  by  Smith,  Woolf  and 
Fleming10.  In  this  method,  the  measured  radiances  l(v)  are  used  to  find  their  equivalents 
at  a single  reference  frequency  v = r through  the  operations 

Vv) 

lr(v) 

i.e.,  we  convert  the  measured  radiances  to  brightness  temperatures  via  the  inverse  Planck 
equation,  and  then  from  the  brightness  temperatures  determine  the  corresponding  radiances 
at  v = r.  (The  reference  frequency  should  be  near  the  median  frequency  of  the  sounding 
channels.)  An  approximate  radiative  transfer  equation  is  obtained  by  writing  (17)  for  the 
reference  frequency  (except  that  the  Tv's  are  unchanged): 

Bo  dTv<p) 

lr(v)  = Br  [T(p0)]  Tv<P0)-/  Br  [T(P,J  SHpf  dx(p)  <l8) 

0 

which  is  linear  in  Br  [T(p)].  (In  fact,  it  is  a linear  Fredholm  equation.)  Once  we  solve  for 
Br  [T(p)]  it  is  a simple  matter  to  obtain  the  temperature  profile. 

Fleming3  has  shown  that  spectral  radiance  errors  resulting  from  the  lineari- 
zation (17)  is  typically  less  than  0.3  percent  over  the  spectral  range  668  cm-1  to  747  cm  . 
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Bv_l  [l(v>], 
Br  [TB(v)]f 


(17) 
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Writing  Equation  18  in  quadrature-perturbation  form  we  obtain 


J 

Alr  (vi}  = Wij  (ABrV  > = I,  2, ...  M,  (19) 

H 


which  in  matrix  notation  is 

Al  = W AB_.  (20) 

~ r ~ r 

Solutions  of  this  equation  are  non-unique,  if  J > M,  and  are  also  unstable  because  of  the  ill- 
conditioned  nature  of  the  weighting  function  matrix  W.  A regularized  solution  of  the  form 

AB  = CAL  (21) 

~ r ~ r 

may  be  obtained  by  requiring  that  Equation  21  represent  a least-squares  best  fit  to  an 

ensemble  of  data,  and  acknowledging  that  the  measurement  vector  components  Al  contain 

ri 

noise  E..  The  least-squares  minimization  procedure  (determination  of  C)  is  completely 
analogous  to  the  derivation  of  Equation  I A,  and  leads  to  the  well-known  matrix  solution 


AB 
~ r 


[’ 


SW'  wsw 


(22) 


where  S is  the  covariance  of  IBr,  and  is  the  M by  M noise  covariance  matrix,  or,  if  the 
channel  noises  are  uncorrelated,  the  diagonal  matrix  of  mean-square  channel  noise  values; 
( )"*"  and  ( )-*  denote  matrix  transpose  and  inverse,  respectively.  Equation  22  is  derived  in 
Appendix  A. 


Not  surprisingly,  there  is  close  similarity  between  Equation  22  and  Equation  I A, 
the  relaxation  formuh  for  the  Fleming  statistical  method.  In  fact,  if  we  truncate  AJr  to  a 
single  component  Alm,  replace  S by  the  matrix  with  components  and  envision  W as  a 
row  vector  of  components  Wm-,  the  two  results  are  identical  expressions  for  AB^;  that  is, 
Equation  22  becomes  a relaxation  formula  suitable  in  a direct  nonlinear  retrieval  method. 


4 


The  minimum  information  equivalent  of  Equation  22  is  obtained  by  assuming 

there  is  no  correlation  between  temperatures  at  different  levels,  and  that  the  reference- 

2 

frequency  Planck  function  variance  is  constant  with  altitude;  i.e.,  S = oB  I,  where  I is  the 

2 r 

identity  matrix.  If  we  also  assume  equal  mean  square  noise  in  the  different  sounding 
channels,  the  Foster  minimum  information  solution  is 

ABr  = WT  [wWT  + oil]"1  Alr;  a = aE2  / aQ  2 (23) 

Equation  23  has  the  same  analogy  to  Equation  3 as  Equation  22  has  to  Equation  14. 

In  the  present  study,  Equations  22  and  23  are  used  to  obtain  highly  over-speci- 
fied retrieval  solutions;  that  is  ABr  is  a vector  of  101  components,  S is  a 101  by  101  matrix, 
W is  6 by  101,  and  the  matrix 

C = SWT  |wSWT  + sj ~ 1 (24) 

is  101  by  6.  The  selected  reference  frequency  is  707  cm  representing  the  fourth  sound- 
ing channel.  The  corresponding  S matrix  is  a subset  of  the  data  S..^  required  by  the 
Fleming  statistical  method,  since  we  also  use  1 01-point  solutions  in  its  evaluation. 

2.3  .ADRATURE  GRID  AND  WEIGHTING  FUNCTIONS 

The  1 01-point  grid  that  we  art  using  to  represent  the  temperature  profile  and  to 

integrate  the  radiative  transfer  equation  corresponds  to  100  atmospheric  pressure  levels 
2/7 

spaced  equally  in  p ' , and  to  the  surface  temperature.  The  weighting  functions  W..  were 

1 1 ' 

determined  using  McClatchey’s  computed  transmittance  profiles.  We  interpolated  these 

to  the  100  atmospheric  levels  and  then  averaged  over  the  nine  different  sets  provided  by 

McClatchey  (representing  nine  different  temperature  profiles)  to  insure  that  the  weighting 

functions  would  not  have  any  noticeable  slope  discontinuities.  In  this  study  the  weights  are 

assumed  invariant  with  atmospheric  temperature. 

Figure  I shows  our  weighting  functions  presented  in  the  conventional  form, 
pressure  versus  -dx^/d  log  p.  (In  our  computations,  the  weights  are  defined  as  transmit- 
tance change  between  successive  gradrature  grid  points;  these  are  given  in  tabular  form  in 
Appendix  B.) 
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Figure  I DMSP  Weighting  Functions  Used  for  Temperature  Retrieval 


To  a very  rough  approximation,  the  direct  nonlinear  methods  obtain  the 
retrieved  temperature  profile  as  a superposition  of  weighting  functions,  when  an  unbiased 
(isothermal)  first  guess  is  used.  Thus,  the  retrieved  profiles  would  generally  not  be 
expected  to  contain  any  structure  sharper  than  the  weighting  functions  themselves,  unless 
the  structure  is  introduced  via  the  initial  guess.  Figure  2 shows  attempts  to  construct 
tropopause-like  features  as  sums  of  two  and  three  weighting  functions.  The  three-function 
sum  gives  an  indication  of  the  achievable  vertical  resolution  at  tropopause  levels. 

Figure  3 shows  the  nine  radiosonde  temperature  profiles  on  which  our  weighting 
functions  are  based.  Except  for  the  Point  Mugu  mid-latitude  profile,  the  profiles  are 
representative  of  the  tropics.  The  very  high  degree  of  correlation  exhibited  by  the  tropical 
profiles  suggests  that  the  statistical  retrieval  methods  may  be  quite  useful,  at  least  within 
a given  climatological  class. 

2 .4  RESULTS  OF  TEMPERATURE  RETRIEVALS  BY  THE  NONST ATISTICAL, 

NONLINEAR  METHODS 

Figures  4 through  7 show  retrievals  obtained  using  synthetic  radiances  and 
Smith's  direct  nonlinear  technique.  In  these  tests  we  selected  the  "actual"  profile  and  the 
initial  guess  from  among  three  different  cases:  the  Point  Mugu  profile  in  Figure  3,  the 
tropical  profile  in  Figure  3 with  the  sharpest  tropopause,  and  an  isothermal  (240  K)  profile. 
We  did  not  include  instrument  noise  in  the  simulation.  Figure  4 shows  that  the  method 
cannot  retrieve  a sharp  tropopause  feature,  starting  from  an  isothermal  guess;  the  error  at 
the  tropopause  is  almost  9 degrees  C.  Similar  errors  result  when  the  actual  profile  is  the 
Point  Mugu  case  (Figure  5),  because  of  the  somewhat  irregular  feature  near  the  100  mb 
level.  Figures  6 and  7 illustrate  the  effects  of  the  initial  guess.  Note  that  the  retrieved 
profile  retains  the  character  of  the  guess  profile  except  at  the  lowest  levels,  where  the 
weighting  functions  are  relatively  large.  Of  course,  the  guesses  used  in  these  last  two 
examples  represent  the  wrong  climatology  and,  therefore,  are  not  fair  tests  of  the  accuracy 
of  the  Smith  method. 

Figures  8,  9 and  10  show  results  obtained  by  the  Fleming  nonlinear  (FLM)*,  the 
alternate  Fleminq  (NEW)*,  and  the  Chow  (CHW)**  algorithms***.  In  each  case  the  actual 

5 In  all  simulations  using  FLM  and  NEW,  the  parameter  a was  set  to  zero. 

**  We  used  K = I in  the  simulations. 

***  Note  that  in  the  error  plots,  positive  errors  appear  to  the  right  of  the  reference  line, 
which  is  a different  convention  than  used  in  Figures  4 through  7. 
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Figure  2 Construction  of  Tropopause-like  Features  Using  Sums  of  Weighting  Functions 
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PRESSURE 


Figure  9 Temperature  Retrieval  Using  the  Alternate  Fleming  Nonlinear 
(NEW)  Method 
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Figure  10  Temperature  Retrieval  Using  the  Chow  Method 
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profile  is  the  tropical  case,  and  the  guess  is  isothermal.  A noticeable  feature  of  the  FLM 
and  NEW  retrieved  profiles  is  their  shapes  in  the  lower  troposphere,  which  closely  resemble 
the  727  cm  * and/or  747  cm  * weighting  functions.  It  appears  that  incorporating  the 
weighting  functions  in  the  relaxation  formulas  can  introduce  a significant  bias  to  the 
solutions  and  resultant  large  errors  between  the  600  mb  level  and  the  surface.  On  the  other 
hand  the  maximum  error  near  the  tropopause  is  smaller  for  the  two  versions  of  the  Fleming 
method  than  for  the  Smith  method.  The  Chow  method  yields  a solution  similar  to  the  one 
obtained  from  the  Smith  method,  with  slightly  greater  error  at  all  levels. 

Figures  8 and  9 show  that  the  retrieved  profile  obtained  by  NEW  is  slightly 
better  than  the  one  obtained  by  FLM,  indicating  that  the  alternate  temperature  averaging 
method  may  yield  a somewhat  less  biased  result. 

Figure  1 1 shows  the  profile  retrieved  by  method  FLM  when  the  actual  profile  is 
the  tropical  case,  and  the  guess  is  the  Point  Mugu  profile  (same  as  Figure  7 except  FLM  is 
used  in  place  of  the  Smith  method  SMH).  Note  that  the  bias  due  to  the  727  cm  * and 
747  cm  * weighting  functions  is  much  less  than  when  an  isothermal  guess  is  used  (Figure  8), 
which  is  apparently  the  result  of  using  a much  more  accurate  guess  in  the  lower 
troposphere  (i.e.,  correct  lapse  rate  but  20-degree  temperature  difference  compared  to  a 
40  to  60-degree  difference  for  the  isothermal  guess).  Note  that  the  initial  guess  imposes  a 
significant  bias  above  the  200  mb  level,  which  is  true  of  all  the  direct  nonlinear  methods. 

Figure  12  illustrates  how  the  accuracy  of  the  retrieved  profile  increases  as  the 
actual  profile  becomes  smoother:  using  one  of  the  smoother  tropical  cases  as  the  actual 
profile,  and  starting  from  an  isothermal  guess,  the  Smith  method  achieves  a retrieval 
accuracy  of  approximately  3 degrees  or  better,  except  at  the  surface  and  tropopause  levels, 
where  the  error  is  about  5 degrees. 

Figures  13  and  14  show  retrieval  results  from  the  Twomey-like  method  TWM 


and  its  alternate  TW2.  Note  that  the  retrieved  profiles  are  similar  to  those  from  the  FLM 
and  NEW  methods,  but  the  errors  are  generally  larger.  TW2  gives  better  retrievals  than 
TWM  below  the  200-mb  level,  but  TWM  does  better  near  the  tropopause. 


TEMPERATURE  (K) 

Figure  12  Temperature  Retrieval  Using  the  Smith  Method 


PRESSURE  (mb) 


Figure  13  Temperature  Retrieval  Using  the  Twomey-Jike  MethtxJ 
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Figure  14  Temperature  Retrieval  Using  the  Alternate  Twomey-like  Method 


The  results  shown  in  Figures  4 and  12  are  indicative  of  the  maximum  accuracy 
achievable  by  the  direct  nonlinear  methods  when  an  unbiased  (isothermal)  first  guess  is 
used,  and  the  actual  profile  is  representative  of  the  tropics.  The  retrieved  profiles  have 
the  general  character  of  a tropical  profile  with  a somewhat  smoothed  tropopause;  hence, 
the  use  of  a guess  profile  which  has  tropical-like  temperatures  and  lapse  rates  but  is  non- 
committal as  to  the  level  or  sharpness  of  the  tropopause  (i.e.,  are  smooth  near  the  tropo- 
pause) would  not  improve  the  retrieval  accuracy.  The  only  way  to  reduce  the  5 to  8°C 
maximum  error  is  to  obtain  more  structural  detail,  which,  it  appears,  is  not  possible  except 
by  introducing  a priori  information.  One  can  utilize  known,  natural  correlations  in  atmos- 
pheric structure,  or  rely  on  an  independently  determined  initial  guess  (e.g.,  a forecast 
profile).  The  first  case  corresponds  to  using  one  of  the  statistical  methods  analyzed  in  the 
present  study;  the  second  represents  the  practice  of  using  the  measured  radiances  merely 
to  effect  corrections  to  a dynamical  forecast  profile. 


Based  on  the  results  given  above,  it  is  our  judgement  that  the  Smith  nonlinear 
method  is  the  "best"  of  the  nonlinear,  nonstatistical  methods  studied.  As  stated  earlier,  the 
inverse  nonlinear  method  of  King"*  was  not  evaluated  in  the  present  study. 


2.5  PROPERTIES  OF  THE  C-MATRICES 


Having  selected  the  Smith  nonlinear  method  as  the  best  representative  of  the 
nonstatistical,  nonlinear  class,  it  remains  to  compare  this  technique  to  the  ones  selected 
for  the  other  three  classes  (Table  3).  These  other  three  classes  use  so-called  C-matrices 
which  together  with  the  measured  radiances  contain  all  of  the  structural  information  avail- 
able. Before  comparing  the  four  methods  or  classes  we  will  examine  the  properties  of  the 
C-matrices. 

The  relaxation  formula  or  matrix  solutions  used  by  the  Fleming  statistical  and 
the  linear  methods  ("full  statistics"  and  "minimum  information")  are 


r(F)  . 

Cji  4li 

(Fleming  Stat.) 

(25) 

- c(fs>ai 

~ r 

(Full  Stat.) 

(26) 

. C(MI)  Al 
~ r 

(Min.  Info.) 

(27) 
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where  the  matrices  C 


(F)  ^(Ml)  anc|  q(FS)  are^  from  Equations  1 4,  23  and  22, 


J 


E 


E 


'28) 


C(MI) 


a 


2.  2 

aE  ^aBr  5 


(29) 


In  Equation  28  the  factor  [ ] is  the  denominator  of  the  expression,  while  in  Equations  29 

and  30  these  factors  are  matrix  inverses.  W is  the  weighting  function  matrix  and  S is  the 

atmospheric  radiance  covariance  matrix  (for  channel  i in  Equation  28,  and  for  the  reference 

or  channel  4 frequency,  v = 707  cm  * , in  Equation  30).  Note  that  the  C-matrices  given  by 

Equations  28  through  30  are  constant  throughout  the  iterative  temperature  retrieval,  since 

W is  assumed  independent  of  temperature.  In  this  study  we  also  assume  that  the  channel 

2 

noises  are  equal  and  uncorrelated,  i.e.,  that  Sp  = aF  I,  where  I is  the  six-by-six  identity 

matrix.  We  used  the  nominal  noise  equivalent  radiance  ac  = 2.5  x 10  W/cm  -sr-cm  = 
2 — | t 

0.25  ergs/s-cm  -sr-cm  for  the  DMSP  temperature  sensing  channels. 

The  tropical  temperature  profile  is  distinctive  in  two  respects:  its  very  sharp 
tropopause  is  difficult  to  reproduce  in  the  retrieved  profiles,  and  the  profile  as  a whole  is 
much  less  variable  than  extratropical  temperature  profiles.  For  this  reason  we  computed 
one  S-matrix  using  data  obtained  at  all  latitudes,  and  another  using  only  tropical  data 
(actually,  the  tropics  was  somewhat  under-represented  in  the  global  data  set,  so  that  the 
former  matrix  may  be  regarded  as  extratropical).  It  is  also  appropriate  to  use  different 
values  for  a (in  Equation  29)  for  the  tropics  and  extratropics.  The  corresponding 
C-matrices,  therefore,  are  dually  defined  and  number  six  in  total. 
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The  average  variance  in  the  global  or  extratropical  reference-frequency  Planck 

radiance  profile,  obtained  by  averaging  the  diagonal  elements  of  the  corresponding 
2 —12  2 —12 

S— matrix,  is  Og  = 1.98  x 10  (W/cm  — sr-cm  ) , which  means  that  the  assumed  noise- 
equivalent  radiance  corresponds  to  a "signal-to-noise"  of  Og  /aE  =*  56  for  the  global  case. 

f I O 

For  the  tropics  the  average  variance  of  the  radiance  profile  is  only  1.42  x 10 
(W/cm^-sr-cm  corresponding  to  a "signal-to-noise"  of  approximately  15.  The 
corresponding  values  of  a are 


2 2 a 

a = Og  /ag  = 3.2  x 10  lExtratropical) 
r 

= 4.4  x 10  ^(Tropics). 


(31) 


We  note  that  the  retrieved  corrections  to  the  guess  temperature  profile  are  constrained  to 
be  smoother  for  the  tropical  case,  but  the  guess  (mean)  profile  has  sharper  structure; 
hence,  the  final  retrieved  profiles  are  not  necessarily  smoother  in  the  tropical  case. 

Equations  26  and  27  show  that  in  the  two  linear  methods  the  solution  vector  AB 
at  each  stage  of  iteration  is  a linear  superposition  of  the  columns  of  or  i.e., 

6 

(4Br>j  * £ Cji  4'r  (vi)' 
i=  I 

(32) 

C.  = C..<FS>orC..(M,) 

J'  Ji  Ji 


Moreover,  since  C is  treated  as  a constant  matrix  in  this  study,  the  AMs  in  Equation  32  may 
be  regarded  as  having  been  summed  over  iteration  number,  and  therefore  AB^  can  represent 
the  total  correction  to  the  initial  guess  Planck  function  profile  13..  In  other  words  the  final 
reference-frequency  Planck  function  profile  is  equal  to  the  guess  profile  plus  a linear 
superposition  of  the  columns  of  C.  Therefore,  a large  set  of  temperature  retrievals 
belonging  to  a given  climatological  class  can  be  economically  stored  or  archived  in  terms  of 
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one  C-matrix,  the  single  initial  guess  profile,  and  six  numbers  per  retrieved  profile,  equal 
to  the  six  Alr(Vj)  summed  over  iteration  number.  In  this  representation  the  summed  Alr 
vectors  are  the  "solutions",  and  the  columns  of  C play  the  role  of  basis  vectors. 

Although  the  influence  of  temperature  in  the  weighting  function  is  ignored  in 
the  present  study,  it  should  probably  be  included  in  the  routine  reduction  of  satellite 
radiance  data.  This  would  require  recomputing  the  W and  C-matrices  in  each  iteration,  and 
would  result  in  retrieved  profiles  that  must  be  represented  by  101  values  (the  number  of 
quadrature  points)  rather  than  six. 

In  the  Fleming  nonlinear  statistical  method,  the  Smith  channel  averaging 

scheme  is  applied  to  the  AB..  obtained  from  Equation  25.  That  is,  the  columns  of  are 

•J  (tr\ 

combined  in  a nonlinear  fashion  to  obtain  the  temperature  profile,  which  means  that  C'  ' is 

not  directly  comparable  to  and  However,  the  following  first-order  analysis  of 

(F) 

the  operations  performed  on  C leads  to  a related  matrix-  K,  whose  columns  are 
approximate  linear  bases  for  the  solution  profiles. 

The  iterate  temperature  profile  T^  obtained  by  evaluating  the  relaxation  Equa- 
tion 25  and  then  averaging  over  channels  is 

6 r l 

Z>ij  Ti  h * ci/F>4li|’ 

T.  = (33) 

J 6 


where  T.(B)  denotes  the  inverse  Planck  function  applied  at  frequency  i,  and  B..  are  the 
vertical  Planck  function  profiles  corresponding  to  the  guess  temperature  profile  (or  pre- 
vious iterate)  ?..  To  first  order,  the  effect  of  a change  AB.j  in  the  argument  of  the  inverse 
Planck  may  be  represented  by 

T.(B..  + AB..)  T.  + AB..  (dBj./dT)-1,  (34) 
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which,  when  applied  to  Equation  23  gives 


t.  = t.  + bl 
j j 


y:  w..  c..( 

Z-/  ij  ji 


(dB../dT)~  Al. 


Z>- 


The  effect  of  a temperature  change  in  the  direct  Planck  function  applied  at  the  reference 
frequency  vf  = 707  cm-'  is,  to  first  order, 

Bfj  (T.  ♦ AT.)  = Brj  (T.)  + AT.  dBrj/dT.  (36) 


Applying  Equation  36  to  Equation  35,  we  obtain  the  reference-frequency  Planck  function 
solution  profile  as 


ri  Ewiicji<F)(<VdT)'l4'i 


B . = B . + 

rJ  r>  dT 


O 

V w.. 

'J 


or  in  matrix  notation 


Br  = ir  + K£l 


By  analytically  evaluating  the  Planck  function  derivatives  in  Equation  37a  we  find  that  the 
elements  of  matrix  K can  be  written 


K [Vf  COSH  (°2  Vf/V  ' ' WiJ  °ji 


IJ  LVi  J cosh  (c-  v /T.)  - I 6 

* ' J 

X,  wii 


where  C2  = 1 .4389  cm°K  is  the  second  radiation  constant. 
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We  see  from  Equation  37b  that  matrix  K plays  the  same  role  as  the  matrices 
C<FS>  or  C(MI>  in  the  linear  methods.  Of  course  the  K-matrix  and  Equation  37b  are  an 
accurate  representation  of  solutions  obtained  by  the  Fleming  statistical  method  only  when 
the  actual  temperature  profile  is  close  to  the  mean  profile  T.  used  for  the  initial  guess. 


It  is  evident  that  plots  of  the  columns  of  C in  the  form  C..  versus  p.  represent  a 
complete  display  of  the  structural  temperature  information  available  in  the  radiance 
measurements,  for  the  two  linear  methods.  Similar  plots  of  the  K-matrices  are  an  approxi- 
mate representation  of  this  information  for  the  Fleming  statistical  method.  The  matrix 
columns  computed  for  the  three  methods  using  both  tropical  and  extratropical  statistics  are 
shown  in  Figures  15  through  20.  For  completeness  we  also  show  the  C-matrices  for  the 
Fleming  statistical  model  (Figures  2 1 and  22),  and  the  mean  reference-frequency  Planck 
radiance  profiles  for  the  tropics  and  extratropics  (Figure  23). 


Figure  15  shows  the  C-matrix  columns  computed  for  the  full-statistics  case 
using  global  (extratropical)  statistics.  A noteworthy  feature,  which  is  common  to  all  of  the 
different  C-matrices,  is  the  channel-to-channel  change  in  vertical  distribution  of  structure 
(information),  which  of  course,  is  due  to  the  influence  of  the  weighting  functions. 
Figure  16  shows  the  minimum  information  C-matrix  for  extratropical  latitudes.  Differ- 
ences between  these  two  cases  are  attributable  to  statistical  temperature  correlations 
between  different  pressure  levels.  Note,  in  particular,  that  the  full  statistics  method  can 
retrieve  more  information  than  the  minimum  information  method  above  the  10-mb  level. 
(Note  the  different  scale  factors  in  Figures  15  and  16.)  The  differences  at  the  lower  levels 
is  essentially  in  small-scale  structure.  Some  of  the  structural  detail  in  the  full  statistics 
matrix  columns  may  be  artificial,  since  we  used  only  100  radiosonde  profiles  in  generating 
the  covariance  matrix.  However,  at  least  some  of  the  structure  exhibited  near  the  300-mb 
level  is  probably  real,  since  there  is  typically  an  abrupt  change  in  temperature  lapse  rate 
near  this  level  at  extratropical  latitudes  (see,  Figure  23,  showing  the  mean  Planck  radiance 
profile).  Correlated  variations  in  the  lapse  rates  just  above  and  below  300  millibars,  or 
variations  in  the  altitude  of  the  tropopause  lower  boundary,  would  be  expected  to  produce 
the  sharp  features  displayed  in  Figure  1 5. 


The  discontinuities  on  the  C-matrix  columns  at  the  lowest  pressure  level  (the 
surface)  arise  from  the  discontinuity  in  the  weighting  functions  at  the  surface  (correspond- 
ing to  the  boundary  term  of  the  radiative  transfer  equation,  which  is  represented  in  the  last 
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Columns  For  The  Full-Statistics  Method  Using  Global  (Extratropical)  Statistics 
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^-Matrix  Columns  For  The  Minimum-Information  Method,  Global  Case  (a=  3.2 


Figure  17  C-Matrix  Columns  For  The  Full-Statistics  Method  Using  Tropical  Statistics 
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Figure  18  C-Matrix  Columns  For  The  Minimum-Information  Method,  Tropical  Case 
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Figure  19  K-Matrix  Columns  For  The  Fleming  Method  Using  Global  Statistics 
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Figure  20  K-Matrix  Columns  For  The  Fleming  Method  Using  Tropical  Statistics 


Figure  21  C-Matrix  Columns  For  The  Fleming  Method  Using  Global  Statistics 


MATRIX  C(F)  (TROPICAL  STATISTICS) 


Columns  For  The  Fleming  Method  Using  Tropical  Statistics 


column  of  W or  last  row  of  W^).  It  may  be  noted  that  the  C-matrix  surface  discontinuities 
would  be  considerably  reduced  in  the  full  statistics  case  if  the  S matrix  exhibited  strong 
correlation  between  the  surface  temperature  and  low-level  air  temperature,  i.e.,  if  SO  00, 
I0l)-«S(I0I,  101).  However,  since  surface  temperature  data  was  not  available  with  the 
radiosonde  data,  we  assumed  no  correlation  between  the  surface  and  any  other  level.  That 
is,  we  set  S(k,  I0l)  = 0;k  = 1,2,...  IOOandS(IOI,  I0I)  = S(I00,  100). 

r Figure  17  applies  to  the  full-statistics  model,  tropical  case,  and  Figure  18  to  the 

minimum  information  model,  tropical  case.  As  expected,  the  C-matrix  columns  are 
smoother  (and  display  smaller  amplitude  modes)  than  in  the  extratropical  cases.  It  is 
interesting  that  the  C-matrix  for  the  full-statistics  tropical  case  (Figure  17),  which  was 
computed  from  175  radiosonde  profiles,  does  not  display  any  sharp  structure  near  the 
tropopause  (100-mb)  level.  Apparently,  it  is  the  gradients  in  the  Channel  I and  Channel  2 
columns  and  the  modal  peak  in  the  Channel  3 column  near  100  mb  that  can  effect  a 
sharpening  or  blunting  of  the  mean  tropical  guess  profile,  shown  in  Figure  23.  The  sharp- 
ness of  the  100-mb  tropopause  in  the  mean  profile  and  the  lack  of  sharp  features  in  the 
C-matrix  columns  near  100-mb  (in  both  the  full-statistics  and  minimum  information  cases) 
shows  that  the  tropopause  will  always  occur  at  the  100-mb  level  in  the  retrieved  profiles. 

Figures  19  and  20  show  the  K-matrices  computed  for  the  Fleming  statistical 
model  using  global  statistics  and  tropical  statistics,  respectively.  Note  that  the  matrix 
columns  are  5 to  10  times  smaller  in  amplitude  than  for  the  other  methods.  This  is  con- 
sistent with  our  observation  that  the  iterative  retrieval  requires  roughly  five  to  ten  times 
more  iterations  before  convergence  is  achieved.  Note  also  that  the  K-matrix  columns  have 
fewer  oscillations  than  the  C-matrix  columns  of  the  other  methods.  This  suggests  (and 
simulated  retrievals  confirm)  that  the  Fleming  statistical  method  is  not  able  to  accurately 
retrieve  profiles  belonging  to  a different  climatological  class  than  the  covariance  matrix 
and  the  mean  profile.  On  the  other  hand  the  method  is  found  to  give  extremely  good 
results  when  appropriate  statistics  are  used. 

2.6  COMPARISON  OF  TEMPERATURE  RETRIEVALS  BY  FOUR  CLASSES  OF 

RETRIEVAL  TECHNIQUE 

In  this  section  we  compare  temperature  retrieval  errors  for  the  Smith  nonlinear 
method,  Fleming  statistical  method,  full-statistics  method  and  minimum  informction 


method.  Comparisons  are  made  separately  for  arctic,  midlatitude  and  tropical  cases. 
When  using  the  last  three  methods  for  either  arctic  or  midlatitude  cases,  we  used  the  mean 
profile,  C-matrices  and  a value  computed  from  global  (extratropical)  data;  the  tropical 
versions  of  these  quantities  were  used  in  tropical  retrievals,  except  where  noted.  The  same 
guess  profiles  (global  or  tropical  mean)  were  used  in  retrievals  by  the  Smith  method.  The 
surface  temperature  is  approximately  equal  to  the  air  temperature  just  above  the  surface 
in  all  cases  (in  both  "actual"  and  guess  profiles). 

2.6. 1 Arctic  Cases 

The  first  arctic  case  (actual  temperature  profile)  is  a 7/22/69  radiosonde  at 
Thule.  Figures  24  through  27  show  the  retrieval  results,  and  Figure  28  compares  the  errors 
for  the  four  methods.  The  Smith  and  Fleming  statistical  methods  give  the  best  overall 
temperature  retrievals;  the  latter  yields  smaller  error  over  the  tropopause  region,  but 
introduces  a small  amount  of  false  structure  at  the  bottom  of  the  tropopause  and  results  in 
a much  larger  error  in  surface  temperature.  Agreement  between  the  Fleming  solution  and 
the  actual  profile  between  0.5  mb  and  10  mb  is  fortuitous. 

A portion  of  the  surface  temperature  error  in  the  Fleming  and  full-statistics 
methods  may  be  attributed  to  the  fact  that  we  could  not  obtain  data  representing  the 
correlation  between  surface  temperature  and  the  temperature  at  different  atmospheric 
levels.  To  examine  the  effects  of  possible  correlations  we  recomputed  the  S-  and 
C-matrices  for  the  F and  FS  methods  using  the  assumption  that  the  surface  temperature  is 
equal  to  the  air  temperature  just  above  the  surface,  (i.e.,  we  set  the  last  row  and  last 
column  of  S equal  to  the  1 00th  row  and  1 00th  column).  This  assumption  eliminates  the 
surface  discontinuity  in  the  full-statistics  solution  and  reduces  it  in  the  Fleming  method. 
The  corresponding  temperature  retrievals  contained  the  errors  shown  by  the  second  and 
third  curves  in  Figure  29  (the  first  and  fourth  curves  are  the  same  as  Figure  28).  Note  that 
the  error  at  the  surface  is  now  roughly  the  average  of  the  previously  calculated  errors  at 
th  surface  and  just  above  the  surface.  A similar  averaging  effect  extends  up  to 
appro  -:mately  800  mb.  Note  also  that  the  Fleming  solution  is  now  slightly  better  than  the 
Smith  solution,  i.e.,  is  the  best  of  the  four  retrievals. 

One  can  generalize  on  the  results  shown  in  Figures  28  and  29:  errors  in  inferred 


surface  temperature  generally  have  compensating  errors  at  some  slightly  higher  altitude 


TEMPERATURE  (1C) 


Figure  25  Simulated  Retrieval  by  the  Fleming  Method.  Actual  Profile 
is  Thule  7/22/69  Radiosonde. 


TEMPERATURE  IK) 


Figure  27  Simulated  Retrieval  by  the  Minimum-Information  Method.  Actual 
Profile  is  Thule  7/22/69  Radiosonde. 


between  the  surface  and  approximately  700  mb.  The  use  of  additional  very  transparent 
channels  to  accurately  determine  the  surface  temperature  would  improve  considerably  the 
accuracy  of  the  retrieved  profiles  below  700  mb,  especially  if  there  are  statistical 
correlations  between  the  surface  temperature  and  low-level  air  temperature.  To  the 
extent  that  data  is  available  for  defining  these  correlations,  it  should  be  used  in  computing 
the  S - and  C-matrices. 

The  second  arctic  case  is  a 11/16/65  radiosonde  at  Heiss  Island  (latitude  81° 
North).  Figures  30  through  33  show  the  retrieval  results,  and  Figure  34  compares  the  errors 
for  the  four  methods.  The  Smith  and  full-statistics  methods  give  the  best  results  for  this 
case.  If  the  temperatures  corresponding  to  levels  100  and  101  are  assumed  perfectly 
correlated,  the  errors  in  the  retrieved  profiles  are  as  shown  in  Figure  35.  Then  the  Fleming 
method  gives  the  best  result  in  our  judgment,  but  it  is  only  marginally  better  than  the 
Smith  or  minimum-information  solutions. 

2.6.2  Midlatitude  Cases 

The  first  midlatitude  case  is  a 2/17/75  radiosonde  at  Point  Mugu.  Figures  36 
through  39  show  the  retrieval  results,  and  Figure  40  compares  the  errors  for  the  four 
methods.  In  this  case  the  full-statistics  method  gives  the  best  result,  with  the  Smith 
solution  running  a close  second.  Figure  41  shows  that  when  the  surface  temperature  is 
assumed  correlated  with  the  above-surface  air  temperature,  the  Fleming  solution  is  best 
and  the  full-statistics  solution  is  second  best. 

The  second  midlatitude  case  is  a 10/22/66  radiosonde  at  West  Geirinish  (latitude 
57°  North).  Figures  42  through  45  show  the  retrieval  results,  and  Figure  46  compares  the 
errors  of  the  four  methods.  The  errors  obtained  assuming  no  surfoce  discontinuity  are 
shown  in  Figure  47.  The  full-statistics  method  gives  the  best  retrieval  in  the  first  case.  It 
also  gives  the  best  retrieval  in  the  second  case,  but  then  it  is  only  slightly  better  than  the 
Fleming  method. 

2.6.3  Tropical  Cases 

The  two  tropical  cases  are  a 2/27/75  radiosonde  at  Kwajalein  and  a 5/8/67 
radiosonde  at  Fort  Sherman.  In  both  cases  the  four  methods  give  nearly  identical  retrie- 
vals, and  those  obtained  by  the  Fleming  and  full-statistics  methods  are  virtually  unaffected 
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Figure  31  Simulated  Retrieval  by  the  Fleming  Method.  Actual  Profile 
is  Heiss  Island  1 1/16/65  Radiosonde. 
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Figure  33  Simulated  Retrieval  by  the  Minimum-Information  Method. 
Actual  Profile  is  Heiss  Island  1 1/16/65  Radiosonde. 


TEMPERATURE  (Is) 


Figure  37  Simulated  Retrievals  by  the  Fleming  Method.  Actual  Profile 
is  Point  Mugu  2/17/75  Radiosonde. 
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Figure  38  Simulated  Retrievals  by  the  Full-Statistics  Method.  Actual 
Protile  is  Point  Mugu  2/17/75  Radiosonde. 


Comparison  of  Retrieval  Errors  for  the  Four  Methods.  Point  Mugu  2/17/75. 
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FOR  CORRELATED  LEVELS  100,  101 


POINT  MUGU 
2/17/75 


Figure  41  Comparison  of  Retrieval  Errors  when  Levels  100,  101  are 

Assumed  Correlated  in  the  F and  FS  Methods.  Point  Mugu  2/ 1 7/75. 
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Figure  45  Simulated  Retrieval  by  the  Minimum-Information  Method. 
Actual  Profile  is  West  Geirinish  10/22/66  Radiosonde. 
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WEST  GEIRINISH 
10/22/66 


Retrieval  Errors  for  the  Four  Methods.  West  Geirinish  10/22/66. 


when  the  alternate  S-matrix  (assumption  of  no  surface  discontinuity)  is  used.  Therefore, 
we  show  only  one  retrieval  (the  one  obtained  by  the  Fleming  method)  for  each  of  the  tro- 
pical cases;  Figure  48  shows  the  Kwajalein  retrieval  and  Figure  49  the  Fort  Sherman  retrie- 
val. The  corresponding  temperature  errors  are  shown  in  Figure  50. 

The  very  different  character  of  the  plots  of  the  K-matrix  columns  of  the 
Fleming  method  (Figure  19  or  20)  and  the  C-matrix  columns  of  the  full-statistics  method 
(Figure  15  or  17)  suggests  that  retrievals  obtained  by  these  two  methods  may  be  affected 
quite  differently  by  the  use  of  inappropriate  statistics.  Figure  51  shows  a retrieval 
obtained  by  the  Fleming  method  for  a tropical  case  when  the  mean  (guess)  profile  and 
C-matrix  are  those  defined  for  extratropical  latitudes.  Figure  52  represents  the  same 
case,  except  that  the  full-statistics  method  is  used  in  the  retrieval.  For  completeness  we 
also  show,  in  Figure  53,  the  results  obtained  when  the  minimum  information  method  is  used. 
Note  that  the  full-statistics  method  yields  a rather  accurate  retrieval  despite  the  use  of 
inappropriate  statistics;  in  fact  the  errors  below  10  mb  have  approximately  the  same 
vertical  distribution  and  are  only  slightly  larger  than  when  correct  statistics  are  used.  The 
retrievals  of  the  other  two  methods  have  unacceptably  large  errors. 

As  remarked  earlier,  the  tropical  mean  profile  and  C-matrices  of  the  linear 
methods  are  such  that  these  methods  will  always  locate  the  tropopause  at  100  mb  in  the 
retrieved  tropical  profiles.  This  is  illustrated  in  Figure  54,  which  shows  the  results 
obtained  using  the  full-statistics  method  when  the  actual  profile  is  a 2/24/75  radiosonde  at 
Barking  Sands.  The  actual  tropopause  is  abnormally  high  (at  approximately  80  mb),  which 
results  in  an  error  of  approximately  8°C  at  this  level  in  the  retrieved  profile. 

2.6.4  Relative  Capabilities  of  the  Retrieval  Methods 

Table  5 summarizes  the  results  described  in  Sections  2.6. 1 to  2.6.3.  The  full- 
statistics  method  provides  the  best  temperature  retrievals  in  the  majority  of  cases  when 
the  S-matrix  is  computed  on  the  assumption  that  the  surface  temperature  is  uncorrelated 
with  the  air  temperature  just  above  the  surface,  although  in  many  cases  it  is  only  mar- 
ginally better  than  the  Smith  method.  When  perfect  correlation  is  assumed,  the  Fleming 
statistical  method  gives  the  best  results  overall,  although  in  two-thirds  of  the  cases  it  is 
equal  to  or  only  marginally  better  than  the  full-statistics  method.  These  results  show  that 
the  full-statistics  method  is  less  sensitive  than  the  Fleming  method  to  the  assumed  air- 
surface  temperature  correlation.  Moreover,  it  is  much  more  tolerant  of  the  use  of 
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Figure  48  Simulated  Retrieval  by  the  Fleming  Method.  Actual  Profile 
is  Kwajelein  2/27/75  Radiosonde. 
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Figure  50  Approximate  Retrieval  Errors  (All  Methods)  For  The 

Tropical  Cases,  Kwaielein  2/27/75  and  Fort  Sherman  5/8/67. 
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Figure  52  Simulated  Retrieval  by  the  Full-Statistics  Method  Using 
Inappropriate  Statistics.  Actual  Profile  is  Kwajelein 
2/27/75  Radiosonde. 
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Figure  53  Simulated  Retrieval  by  the  Minimum-Information  Method  Using 
Inappropriate  Statistics.  Actual  Profile  is  Kwajelein  2/27/75 
Radiosonde. 
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Figure  54  Simulated  Retrieval  by  the  Full-Statistics  Method.  Actual 
Profile  is  Barking  Sands  2/24/75  Radiosonde. 
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Table  5 METHODS  GIVING  BEST  TEMPERATURE  RETRIEVALS  FOR  ARCTIC, 
MIDLATITUDE  AND  TROPICAL  CASES 


Best* 

Best** 

Max 

Error*** 

TO 

Comments 

ARCTIC 

Thule 

SMH  & F 

F 

2.5 

Heiss  Island 

SMH  & FS 

F(sFS) 

4 

MID- 

LATITUDE 

Pt.  Mugu 

FS 

F 

3 

West  Geirinish 

FS 

FS(=F) 

2.5 

TROPIC 

Kwajalein 

4.5 

SMH,  F,  FS  & Ml 

Fort  Sherman 

— 

— 

2.5 

give  equal  results 

* S-Matrix  assumes  no  correlation  between  Levels  1 00,  101. 
**  S-Matrix  assumes  equal  temperature  at  Levels  100,  101. 
***  In  solution  "Best**"  at  Levels  below  70  mb. 


inappropriate  statistics  in  defining  the  other  S-matrix  elements,  which  suggests  it  will  give 
better  results  in  the  case  of  anomalous  temperature  profiles.  On  the  basis  of  these 
considerations,  we  believe  the  full-statistics  method  is  the  best  of  the  nine  different  linear 
and  direct  nonlinear  techniques  studied,  at  least  for  the  retrieval  of  temperature  structure. 

Maximum  errors  below  the  70  mb  level  at  extratropical  latitudes  are  approxi- 
mately 3 to  4 degrees  C when  noise  in  the  data  is  negligible.  For  the  tropics  the  maximum 
error  is  between  4 and  5 degrees,  except  when  the  tropopause  is  unusually  high  or  low,  when 
much  larger  errors  will  occur. 


82 


7809-4 


SECTION  3 

STUDY  OF  WATER  VAPOR  RETRIEVAL 
3. 1 SELECTION  AND  DEFINITION  OF  ALGORITHMS 

The  radiative  transfer  equation,  usually  written  in  the  form 

t(v,  ps) 

l(v)  = B(v,  Ts)  t(v,  ps)  — j B (v,  p)  d t(v,  p),  (39) 

I 

can  be  integrated  by  parts,  which  gives 


l(v)  = J^B  (v,  T#)  - B (v,  p$)j  t(v,  ps) 


(40) 

B(v,ps) 

+ B (v,  pf)  + / t(v,  p)  dB  (v,  p), 

B (v,  pf) 

where  I (v)  is  the  upwelling  spectral  radiance,  t(v,  p)  is  the  transmittance  of  the  atmos- 
phere above  level  p,  and  B (v,  p)  is  Planck's  function  evaluated  at  temperature  T (p); 
Subscripts  s and  t refer  to  the  bottom  and  top  of  the  atmosphere,  respectively;  however,  in 
the  first  term  of  Equation  40,  B (v,  T ) denotes  the  emission  of  the  surface  at  temperature 

T , whereas  B (v,  p ) is  Planck's  function  evaluated  at  the  air  temperature  T (p  ) just  above 
s s s 

the  surface.  It  is  often  assumed  that  T = T (p  ),  which  would  eliminate  the  first  term. 

5 5 


Replace  I (v)  in  Equation  40  by  Al  (v)  and  t(v,  p)  by  A(ln  u*)  dx/dln  u*,  where 
Al  (v)  is  the  perturbation  in  radiance  caused  by  a perturbation  Ain  u*(p),  and  u*(p)  is  defined 
as  the  reduced  water  vapor  optical  thickness  above  level  p: 


This  quantity  is  proportional  to  u(p)  7 (p),  the  actual  water  thickness  times  the  Curtis- 
Godson  mean  line  halfwidth.  The  choice  of  u#(p)  or  In  u*(p)  as  the  variable  defining  the 
atmospheric  water  distribution  is  based  on  the  fact  that  the  DMSP  water  vapor 
transmittances  are  in  the  strong-line  absorption  region  at  nearly  all  pressure  levels  (see 
Section  3.2). 


Writing  Equation  40  in  a quadrature  form  and  in  terms  of  the  perturbations 
Al  (v)  and  A In  u*,  we  obtain 


J 

Al.  Pjj  Ainu*;  i = I,  2, ...  6 (42) 

j=l 


where 


P. 


dTi  I r 1 

B - Bj  | 
dlnu*  I 


I 


P - - 
*,J  2 din 


[Bi,j  + I Bi,j-l]  » j = I • (43) 


’i.J  = 1 “d]^5  [^I.J  + I 8i,J  Bi,J  - |] 
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where  J is  the  number  of  quadrature  points  in  the  atmosphere,  and  B.  is  the  surface 

i,  J + I 

radiance.  As  in  the  analysis  of  temperature  retrieval  techniques,  we  chose  quadrature 

2/7 

points  corresponding  to  equal  increments  of  p . The  mesh  extends  from  the  100  mb  level 
to  the  surface,  and  over  this  region  matches  the  mesh  used  for  temperature  retrievals. 
This  results  in  a total  of  J = 53  quadrature  points.  From  Equations  42  and  43,  we  see  that 
the  vertical  weighting  function  on  Ain  u*  (p^)  is  (dr /din  u*)  (dB/dp^). 

The  different  water  vapor  retrieval  techniques  to  be  examined  will  solve  for  the 
profile  Ain  uf,  representing  the  difference  between  the  actual  In  u*  profile  and  an  initial 
guess  profile;  Al.  are  differences  between  the  measured  radiances  and  those  computed  for 
the  guess  profile.  In  subsequent  iterations  the  solution  profile  replaces  the  initial  guess  or 
previous  solution.  The  solution  has  converged  when  (and  if)  the  Al.'s  become  less  than  the 
rms  noise  level.  The  profile  of  integrated  water,  u(p),  can  be  obtained  from  u*(p)  using  the 
inverse  of  Equation  41: 


(44) 


The  water  vapor  algorithms  selected  for  study  are  the  same  four  that  were 
studied  with  respect  to  temperature  retrieval  in  Section  2.6.  However,  only  the  two  linear 
methods  were  actually  evaluated  (see  Section  l.l). 

3.1.1  Linear  Inverse  Methods 


In  matrix  notation  Equation  42  is 


Al  = P . (Ainu*) 


and  its  regularized  inverse  is 


Ain  u*  = C A I 


(45) 


(46) 
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If  the  Foster-Twomey  regularization  is  used, 


C = C(M,)  = PT  £pPT  + al  J , 


2,  2 

a = aE  /aln  u*’ 


and  Equation  46  is  the  minimum  information  solution.  Here  I is  the  6-by-6  identity  matrix, 

2 

a£  is  noise  equivalent  radiance  and  <T|n  ^ is  the  expected  variance  of  u*/<u*>  averaged 
over  all  pressure  levels. 


If  statistical  regularization  is  used 


C = C(FS)  = SPT. 


+ SE  ] 


where  S is  the  covariance  matrix  of  u*/<u*>  and  SP  is  the  6-b y-6  matrix  of  noise  covari- 

2 —ft  2 — I 

ances.  We  will  assume  that  |,  with  = 2 x 10  W/cm  -sr-cm  . These  solu- 

tions are  analogous  to  those  used  in  the  linear  methods  for  temperature  retrieval:  Ain  u* 
plays  the  role  of  ABf  and  P the  role  of  W. 

3. 1 .2  Smith  Iterative  Technique 

The  iterative  method  described  by  Smith*  and  also  by  Smith  and  Howell**, 
starts  with  the  assumption  that  the  altitude  variation  of  Ainu*  in  Equation  42  will  be 
relatively  weak  compared  to  that  of  the  weighting  functions  (rows  of  P).  If  Ainu*  is 
assumed  constant  within  the  layers  sensed  by  a particular  channel,  we  can  write  (from 
Equation  42), 


Ain  uj* 


vE  pij 

j=i 


7809-4 


This  is  a relaxation  formula  for  revising  In  uf,  which  yields  six  independent  new  estimates 
of  the  profile  u*  Smith  calculates  a single  profile  as  a weighted  sum  of  the  six  estimates, 

where  the  weighting  function  is  P..: 

ij 


Ain  u* 
J 


6 

E 

"pn  4li  1 

J 

i=  1 

L H J 

(51) 


Note  that  Equation  51  is  a linear  transformation;  in  matrix  notation 

Ainu*  = CA^  (52) 


with 


However,  the  Smith  method  (and  the  "linear"  methods  and  the  Fleming  statistical  methods) 
are  linear  only  through  the  first  solution  of  Equation  52.  In  subsequent  iterations  P and  C 
must  be  recalculated,  since  P..  depends  on  the  current  u*  through  dx.(u^/d  In  u*. 

3. 1 .3  Fleming  Statistical  Method 


As  in  the  linear  and  Smith  iterative  methods  the  solution  has  the  general  form 
of  Equation  52.  The  C-matrix  can  be  obtained  by  analogy  to  the  equations  used  for 


temperature  retrieval:  substitute  Equation  28  into  25  and  then  average  over  the  channels, 
as  in  the  Smith  method.  The  result  is 


(F) 


>..  T,  s.  p. 

Ij  immJ  jn  in 
n=  I 


• • — V— ' • • — 

Jl  Jl 

r~l~7 7” 

E K E 

L i=l 

L J= 1 \ n=l 

(54) 


where  S.  , the  covariance  of  u*/<u*>,  replaces  S..  in  Equation  28,  and  P„  replaces  W... 
Jn  ijn  ij  ^ ij 

3.2  TRANSMISSION  FUNCTIONS  FOR  THE  DMSP  WATER  VAPOR  CHANNELS 


Computed  water  vapor  spectral  transmittance  profiles  t(v,  p)  were  provided  by 
Dr.  McClatchey  of  AFGL  for  218  selected  radiosondes.  Also  provided  were  the  spectral 
responses  of  the  six  DMSP  channels  listed  in  Table  6.  From  these  data  we  computed  218 
transmission  profiles  for  each  of  the  six  channels.  The  results  were  plotted  in  various 
forms  in  order  to  obtain  empirical  transmittance  functions  for  the  six  channels. 


Table  6 DMSP  HzO  CHANNELS 


Channel  Number 
(Used  in  this  Study) 


Nominal  Center 
Frequency 

(cm  ’) 


1 

2 

3 

4 

5 

6 


353 

347 

397 

408 

420 

441 
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Figure  55,  from  Zachor  , shows  the  general  behavior  of  absorptance  versus 
optical  thickness  u and  pressure  p for  a homogeneous  path.  This  behavior  is  typical  of 
laboratory  data  and  the  predictions  of  all  band  models.  Figure  56  shows  computed 
absorptances  for  Channel  4 plotted  in  the  form  of  Figure  55.  Note  that  six  distinct  values 
of  u are  represented,  and  that  the  plotted  points  tend  to  fall  on  a single  curve.  The  plot 
shows  that  most  of  the  data  is  in  the  strong-line  region  to  within  a few  percent  in 
absorptance.  If  most  of  the  Channel  4 transmittances  are  in  the  strong-line  region,  then  so 
will  be  the  transmittances  of  the  other  channels,  since  Channel  4 is  the  most  transparent. 

Figure  57  shows  all  of  the  computed  transmittances  for  Channel  4 plotted  as  log 
(I  - t)  versus  log  u*.  Again,  the  departure  of  the  plotted  values  from  a single  curve  is 
quite  small.  Some  of  the  scatter  is  due  undoubtedly  to  the  effect  of  temperature  in  the 
transmittance.  The  scatter  at  very  low  u*  is  due  to  truncation  of  the  t's  and  u's  to  four 
digits. 

Figure  57  and  similar  plots  for  the  other  five  channels  were  smoothed  and  used 
to  obtain  approximate  empirical  transmission  functions  t(u*)  in  table  form.  We  then  used 
the  empirical  functions  to  obtain  plots  of  dx/d  In  u*  versus  u*,  which  is  a factor  in  the 
weighting  function  referred  to  in  Section  3.1.  These  plots,  however,  exhibited  ragged  fine 
structure  (Figure  58)  which  is  not  physically  plausible  and  could  tend  to  produce  false 
structure  in  the  retrieved  humidity  profiles.  We  then  decided  to  fit  empirical  functions  to 
the  data  in  order  to  obtain  a smooth  representation  of  the  weighting  functions.  For  the 
empirical  transmittance  functions,  we  adopted  the  generalized  strong-line  form 


— In  x = 


KjX 


«2x 


KnX 


x = (u*)K 


(55) 


and  determined  the  coefficients  K.  and  exponent  k with  the  assistance  of  least-squares 
curve  fitting  codes. 
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log  (up) 


Figure  55  A Plot  of  the  log  of  absorption  vs  log  (up)  for  a typical 
laboratory  run  or  absorption  band  model.  The  ordinate 
scale  is  greater  than  the  scale  of  the  abscissa  to  emphasize 
the  strong-line  and  Beer's  law  regions.  The  straight 
portion  of  the  strong-line  envelope  represents  the 
"square-root"  approximation,  A = const.(up)” 
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It  was  found  that  the  empirical  analytical  functions  did  not  give  acceptable  fits 
to  the  transmittance  data  for  high  transmittances.  This  probably  should  have  been 
anticipated,  since  according  to  Equation  55  the  derivatives  d In  t /du*and  dx/du* approach 
zero  as  u*-»0  if  k<  I,  and  approach  infinity  as  u*->0  if  k<  I.  The  corresponding 
weighting  functions  were  well-behaved  and  very  smooth,  but  were  generally  too  low  at 
small  u*  and  were  too  sharply  peaked,  based  on  the  original,  ragged  weighting  functions. 

In  view  of  the  above  difficulties  it  was  decided  to  retain  the  empirical 
transmittance  functions  in  tabular  form,  but  apply  a numerical  smoothing  procedure.  This 
consisted  of  convolving  x^  versus  In  u*^  (the  original  tabular  values)  with  a unit  area 
triangle  of  base  width  Ain  u*  = 0.35.  From  the  resulting  values  we  constructed  a 
corresponding  table  of  derivatives  dx^/d  In  u*.  These,  the  final  weighting  function  values, 
are  plotted  in  Figure  59.  The  smoothed  transmittances  are  shown  in  Figure  60. 

3.2. 1 Transmission  Functions  for  Alternate  Channels 

Alternate,  more  transparent,  channels  would  be  desirable  for  retrieving  low- 
level  moisture,  as  pointed  out  in  Section  4.2.  A transmission  function  for  Channel  E7  at 
535  cm  * could  probably  be  obtained  in  the  same  manner  as  the  transmission  functions 
described  above,  since  the  line-by-line  spectra  provided  by  McClatchey  includes  most  of 
this  band. 

13  -I 

As  suggested  by  McClatchey  , window  channels  near  898  and  797  cm  may  be 

useful  for  inferring  both  total  water  and  surface  temperature,  as  well  as  the  low-altitude 
moisture  distribution  when  these  channels  are  used  in  combination  with  the  existing 
channels.  For  these  channels  McClatchey  suggested  using  the  transmission  function  pro- 
posed by  Goody  * \ 


TRANSMITTANCE  DERIVATIVES 
(FINAL,  SMOOTHED  TABLE) 


1 353  cm' 

2 347 

3 397 

4 408 

5 420 

6 441 


where  the  line  intensity  and  root-intensity-halfwidth  sums  are  over  the  lines  in  the  interval 

Kt 

Av.  McClatchey  has  computed  Y'S.  and  Y'(S'r).  as  a function  of  temperature  for  various 

—I  1 —I 

spectral  intervals,  including  20  cm  bands  at  795  and  900  cm  . Equation  56  does  not 
include  the  contribution  of  the  water  vapor  continuum. 

Equation  56  applies  to  homogeneous  paths.  Following  Godson*"*  one  can  show 
that  an  appropriate  generalization  of  Equation  56  for  an  inhomogeneous  path  is 


t = exp 


-/  £ s,  du 


L.  [j 

S.)  du 

2 z 

* ) 

[£<sL 

i 

du 

(57) 


This  transmission  function  together  with  the  LOWTRAN  4 continuum  model  was  used  for 
computations  involving  the  795  and  900  cm  * channels. 

The  alternate  window  channels  at  795  and  900  cm  * and  the  E7  band  were  not 
used  in  the  formal  retrievals  described  in  Section  3.4.  The  effectiveness  of  the  window 
channels  in  determining  total  water  and  surface  temperature  is  discussed  in  Section  3.5. 

3.3  C-MATRICES  FOR  HjO  RETRIEVAL 

Figures  6 1 and  62  show  C-matrices  computed  for  the  minimum  information 

2 

method.  The  first  corresponds  to  a wet  tropical  profile  with  4.4  gm/cm  total  water,  and 

2 

the  second  to  a dry  arctic  case  with  0.06  gm/cm  total  water.  These  matrices,  which  are 
typically  iterated  several  times,  do  not  characterize  the  retrieval  solutions  in  the  same  way 
as  the  constant  C-matrices  for  temperature  retrieval.  However,  they  indicate  the  pressure 
ranges  over  which  humidity  information  is  retrievable  under  very  wet  and  very  dry 
conditions.  For  the  wet  tropical  case  the  range  is  approximately  200  mb  to  700  mb;  for  the 
dry  arctic  case  it  is  approximately  300  mb  to  the  surface. 
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RESULTS  OF  H20  RETRIEVALS 


3.4. 1 Nature  of  Solutions  for  Tropical  Cases  (Minimum-Information  Solution) 

Figure  63  shows  a simultaneous  water  vapor  and  temperature  retrieval  for  a 
tropical  atmosphere  obtained  by  the  minimum-information  method.  Parts  (a)  through  (d) 
show  retrieved  relative  humidity,  retrieved  log  u(p),  the  corresponding  u(p),  and  retrieved 
temperature,  respectively.  Part  (b)  of  the  figure  best  illustrates  the  nature  of  the  water 
vapor  retrieval  solution. 

The  guess  profile  for  In  u is  based  on  a guess  of  60  percent  relative-  humidity  at 

all  altitudes  (Figure  63  a)  and  the  retrieved  temperature  profile  (Figure  63  d).  The 

2 

resulting  total  precipitable  water  in  the  guess  profile  is  3.7  gm/cm  compared  to 
2 

4.4  gm/cm  in  the  actual  profile.  Note  that  the  retrieved  In  u profile  does  not  depart 
significantly  from  the  guess  above  the  200  mb  level  or  below  700  mb,  due  to  the  small 
amplitude  of  the  weighting  functions  (and  C-matrix  columns)  in  these  regions.  Between 
200  and  700  mb  the  In  u solution  oscillates  about  the  true  profile,  resulting  in  positive  and 
negative  compensating  errors;  the  error  follows  on  oscillatory  pattern  not  unlike  the 
original  C-matrix  column  for  Channel  I (Figure  61).  The  solution  is  the  result  of  three 
iterations  (which  produced  convergence  between  the  measured  and  calculated  radiances). 

Because  the  solution  for  In  u contains  an  oscillatory  component,  the  correspond- 
ing solution  for  relative  humidity,  which  is  proportional  to  the  derivative  of  u,  also  has  an 
oscillatory  component.  In  fact,  the  oscillations  in  relative  humidity  are  relatively  much 
larger,  as  shown  in  part  (a)  of  Figure  63.  Of  course,  our  purpose  is  to  retrieve  the  In  u 
profile.  Had  we  sought  a solution  for  relative  humidity  it  is  possible  that  better  behaved 
solutions  would  have  been  obtained  at  the  expense  of  degraded  In  u solutions. 

Figures  63  a and  c show  that  no  radiation  is  received  from  levels  below  approxi-  , 

mately  700  mb  because  of  the  large  amount  of  total  water,  and  that  approximately 

75  percent  of  the  water  is  below  700  mb.  It  is  therefore  not  surprising  that  the  retrieved 

2 

profile  contains  the  same  total  water  as  the  guess  profile  (~  3.7  gm/cm  ). 


RELATIVE  HUMIDITY  (percent) 


Figure  63a  Water  Vapor  Retrieval  For  A Tropical  Atmosphere  By  The 

Minimum-Information  Method.  (Retrieved  Relative  Humidity). 


PRESSURE  (mb) 


PRESSURE  $nb) 


PRESSURE  (mb) 


Figure  64  shows  the  relative  humidity  retrieved  for  a similar  case.  Here  the 

guess  profile  is  a linear  approximation  (in  relative  humidity  versus  log  p)  to  the  actual 

2 

profile.  The  total  retrieved  water  is  4.5  gm/cm  , approximately  the  same  as  for  the  guess 
and  actual  profiles.  The  retrieved  humidity  profile,  while  much  smoother,  bears  only  a 
slight  resemblance  to  the  actual  profile. 

3.4.2  Channel  Redundancies 

Figure  59,  which  shows  the  weighting  function  factor  -dt/d  In  u*  for  the  six 
channels,  suggests  that  Channel  2 is  somewhat  redundant  with  Channel  3,  and  that  Chan- 
nels 4,  5 and  6 are  redundant  with  each  other.  Since  redundancy  among  channels  tends  to 
promote  instability,  it  is  logical  to  ask  whether  some  of  the  oscillations  in  the  retrieved 
moisture  profiles  are  a result  of  the  redundancies.  We  therefore  repeated  one  of  the 
tropical  retrievals  described  above  using  only  the  three  most  independent  channels:  from 
Figure  59  these  are  Channels  I,  2 and  4,  at  353  cm-*,  347  cm-*  and  408  cm-*,  respec- 
tively. The  conditions  were  the  same  as  those  represented  in  Figure  63. 

v 

The  moisture  profile  retrieved  using  three  channels  differed  only  slightly  from 
that  obtained  using  the  original  six  channels.  In  fact  plotted  in  the  form  In  u versus  In  p, 
the  result  would  be  indistinguishable  from  the  "retrieved"  curve  in  Figure  63  b.  The 
corresponding  relative  humidity  profiles  for  the  3-  and  6-channel  retrievals  are  shown  in 
Figure  65.  Note  that  the  3-channel  result  is  slightly  smoother  than  the  6-channel  result; 
however,  the  differences  in  relative  humidity  values  are  very  small. 

We  also  performed  3-  and  6-channel  retrievals  for  a moderately  dry  case, 
represented  by  a Fort  Churchill  (latitude  59°  N)  radiosonde  with  0.78  gm/crn^  total  preci- 
titable  water.  The  3-  and  6-channel  results  for  both  In  u and  relative  humidity  were  the 
same  within  a few  percent. 

The  above  comparisons  show  that  virtually  al’  the  information  available  in  the 
DMSP  measurements  is  contained  in  Channels  I,  2 and  The  retrievals  presented  in  the 
following  sections  used  only  these  three  channels,  except  where  noted. 
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3.4.3 


Solutions  Obtained  by  the  Full-Statistics  Method 


The  moisture  retrievals  reported  above  were  all  obtained  using  the  minimum 
information  method.  We  performed  additional  retrievals  (for  other  cases)  by  both  the 
minimum- information  and  full-statistics  methods,  and  found  that  the  latter  method 
generally  gives  better  results.  The  retrievals  presented  below  will  be  based  on  the  full- 
statistics  method  except  where  noted.  (This  method  is  used  for  both  temperature  and 
moisture  retrieval). 

Figure  66  shows  the  results  of  a water  vapor  retrieval  for  a moderately  dry, 
2 

high-latitude  case  (0.78  gm/cm  of  water  in  the  actual  profile).  The  guess  is  a mean 
moisture  profile  calculated  from  non-tropical  radiosonde  data.  Instead  of  retrieving  the 
temperature  profile  and  using  it  in  the  h^O  retrieval,  we  used  the  "actual"  temperature 
profile.  These  rather  favorable  conditions  — dry  atmosphere,  reasonable  guess  for  In  u(p), 
no  temperature  retrieval  errors  — would  be  expected  to  result  in  solutions  that  represent 
an  upper  limit  to  the  capabilities  of  the  DMSP  H-package. 

The  retrieved  In  u profile  is  a fairly  good  representation  of  the  actual 

2 

(Figure  66  a).  The  total  precipitable  water  in  the  retrieved  profile  is  1.04  gm/cm  , 
compared  to  0.78  and  0.40  for  the  actual  and  guess  profiles.  The  relative  humidity  profile 
corresponding  to  the  retrieved  In  u profile  (Figure  66  b)  shows  considerable  false  structure, 
and  contains  large  errors. 

Some  of  the  errors  in  relative  humidity  arise  from  the  oscillatory  components  of 
the  retrieved  In  u profile.  This  suggests  that  it  may  be  possible  to  reduce  the  humidity 
errors  by  smoothing  , after  the  fact,  the  In  u solution.  To  test  this  idea  we  replaced  the 
retrieved  In  u versus  In  p profile  of  Figure  66  a by  a least-squares  straight-line  approxima- 
tion. The  corresponding,  smoothed,  relative  humidity  profile  is  shown,  together  with  the 
actual  profile  in  Figure  66  c;  the  former  is  labeled  "Retrieved  (LSQ)".  Also  shown  is  the 
error  in  this  relative  humidity  profile  and  the  corresponding  error  in  integrated  water  u(p); 
these  are  compared  to  the  errors  that  existed  prior  to  smoothing  the  retrieved  In  u(p). 

It  is  evident  that  the  least-squares  smoothing  process  improved  the  inferred 
relative  humidity  distribution,  but  degraded  the  In  u(p)  or  u(p)  profile;  the  third  set  of 
curves  of  Figure  66  c shows  that  smoothing  substantially  increased  the  average  error  in 
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retrieved  u(p).  It  is  noteworthy  that  the  original  In  u(p)  solution  satisfies  the  radiance 
convergence  criterion  but  the  smoothed  version  does  not.  This  example  points  out  that  the 
capabilities  of  a particular  retrieval  algorithm  may  be  quite  different  with  respect  to  point 
variables  (such  as  relative  humidity)  and  integrated  quantities,  like  u(p).  As  stated  earlier 
the  algorithms  considered  in  this  study  were  developed  specifically  for  retrieving  In  u(p). 

Although  the  least-squares  smoothing  improved  the  humidity  retrieval  in  this 
particular  example,  it  is  not  a recommended  procedure.  If  the  actual  In  u versus  In  p 
profile  deviates  significantly  from  a straight  line,  then  forcing  this  behavior  on  the 
retrieved  In  u(p)  can  result  in  a completely  erroneous  distribution  for  relative  humidity.  A 
better  general  procedure,  short  of  using  a completely  different  algorithm  for  relative 
humidity,  would  be  to  simply  increase  the  estimate  of  the  channel  noise  in  Equation  49 
(or  increase  the  parameter  a in  the  minimum  information  method)  when  relative  humidity 
will  be  derived  from  the  u(p)  solution. 

For  the  case  represented  in  Figure  66  we  also  performed  a simultaneous 
temperature-humidity  retrieval;  i.e.,  we  included  the  effects  of  temperature  retrieval 
errors,  which  were  between  +I°C  below  the  400  mb  level  and  between  +3°C  above  400  mb. 
These  errors  produced  small  changes  in  retrieved  u(p)  and  relative  humidity,  but  no 
significant  degradation  in  the  overall  quality  of  the  solutions. 

Figure  67  shows  the  results  of  a tropical  water  vapor  retrieval.  The  actual 

2 2 2 
profile  has  4.4  gm/cm  water,  the  guess  about  I gm/cm  , and  the  retrieved  3.3  gm/cm  . 

The  retrieved  temperature  profile  (not  shown)  has  errors  between  +3°C.  As  Figure  67  a 

shows,  the  retrieved  In  u(p)  profile  is  quite  a good  representation  of  the  actual  profile.  This 

may  be  fortuitous,  i.e.,  may  result  from  the  fact  that  the  guess  and  actual  profiles  are 

similar  in  shape  (although  in  magnitude  they  are  quite  different).  Note  that  the  solution  is 

not  constrained  to  follow  the  guess  below  700  mb  even  though  all  three  channels  are  opaque 

to  radiation  from  below  700  mb;  this  results  from  using  statistical  regularization. 

Figure  67  b shows  the  retrieved  relative  humidity,  which,  between  200  and  700  mb,  is  an 

excellent  representation  of  the  actual  relative  humidity. 

Figure  68  displays  the  least  encouraging  of  the  humidity  retrievals  obtained  in 

this  study.  The  actual  profile  is  a West  Geirinish  (64°  latitude)  radiosonde  with  a rather 

2 

large  amount  of  precipi table  water  (3.6  gm/cm  ).  The  guess  is  the  same  profile  used  as  a 
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200  H- 


Figure  67b  Water  Vapor  Retrieval  For  A Tropical  Case,  (Full-Statistics 
Method).  Retrieved  Relative  Humidity  Profile. 
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guess  in  the  previously  described  retrievals  based  on  the  full-statistics  method:  it 

2 

represents  a mean  for  tropical  latitudes  but  is  scaled  to  I gm/cm  total  precipitable  h^O. 
The  covariance  matrix  S also  represents  tropical  conditions.  The  retrieved  profile  contains 
8.3  gm/cm^  of  precipitable  water. 

The  retrieved  In  u(p)  profile  is  actually  a reasonable  low-resolution  approxima- 
tion to  the  guess  profile.  However,  it  is  much  too  smooth,  which  is  a property  imposed  by 

the  guess  profile  and  the  covariances  used  in  regularizing  the  solution.  The  smoothness  of 

2 

the  solution  results  in  the  rather  large  error  of  4.7  gm/cm  in  total  precipitable  water.  The 
errors  in  u(p)  are  equally  large  on  a percentage  basis  at  the  500  mb  level.  (Note  that  if  we 
had  used  the  minimum  information  method  the  total  water  would  be  determined  by  the 
guess,  and  the  error  would  still  be  large.) 

We  also  computed  a mean  In  u profile  and  covariance  matrix  S for  extratropical 
latitudes.  The  data  base  is  a set  of  1200  radiosondes,  provided  by  AFGL  — the  same  data 
used  for  the  calculation  of  the  tropical  mean  and  covariance.  The  two  means  were  found  to 
have  virtually  the  same  shape;  i.e.,  <u>  for  the  tropics  is  approximately  a factor  of  two 
higher  than  the  extratropical  <u>  at  all  pressure  levels.  It  was  found  also  that  the  water 
vapor  retrievals  obtained  by  the  full-statistics  method  depend  only  slightly  on  which  mean- 
covariance  set  is  used.  Similarly,  it  was  found  that  the  solutions  are  unaffected  when  the 
guess  u(p)  profile  is  scaled  by  a constant  factor.  Thus,  the  errors  in  the  retrieved  profile 
shown  in  Figure  68  cannot  be  attributed  to  using  statistics  for  the  wrong  latitudes,  or  to 
applying  too  small  a scale  factor  to  the  guess  profile. 

3.4.4  Effects  of  Temperature  Retrieval  Errors 

Errors  in  the  retrieved  temperature  profile  will  naturally  contribute  to  errors  in 
the  retrieved  In  u profile.  This  error  propagation  effect  was  not  studied  in  detail;  a few 
sample  computations  were  performed  to  estimate  the  magnitude  of  the  effect  and 
determine  general  trends. 

Figure  69  shows  the  results  of  one  such  calculation.  We  first  retrieved  the 
temperature  profile,  which  contained  errors  between  +2.5°C,  and  a surface  temperature 
error  of  -l.0°C.  We  then  performed  a humidity  retrieval  using  in  one  case  a correct  guess 
for  the  vertical  distribution  of  relative  humidity,  and  in  another  case  a correct  guess  for 


the  precipitable  water  profile  u(p).  When  one  of  these  is  correct  the  other  will  be  incorrect 

due  to  the  temperature  profile  errors.  The  second  two  sets  of  curves  in  Figure  69  shows 

the  resulting  errors  in  the  retrieved  u(p)  profiles  and  the  corresponding  errors  in  relative 

2 

humidity.  Note  that  the  magnitude  of  the  error  in  u(p)  is  of  the  order  of  0.1  gm/cm  , which 
is  small  compared  to  the  retrieval  errors  noted  in  earlier  sections.  This  particular 
calculation  used  the  minimum  information  method  and  a 6-channel  retrieval. 

Sometimes  large  errors  will  occur  in  the  retrieved  surface  temperature  (e.g., 

because  the  actual  surface  temperature  and  -air  temperature  above  the  surface  represent  a 

discontinuity  much  larger  than  that  of  the  guess  profile  or  those  corresponding  to  the  the 

C-matrix  columns).  Then  there  will  tend  to  be  large  compensating  errors  at  low  levels  in 

the  retrieved  temperature  profile,  and  correspondingly  large  errors  in  the  retrieved  In  u 

profile  at  low  levels.  This  effect  was  encountered  in  a 6-channel  arctic  retrieval,  using  the 

minimum  information  method.  The  retrieved  surface  temperature  was  I0°C  too  low,  which 

2 2 

resulted  in  a retrieved  total  precipitable  water  of  0.25  gm/cm  , compared  to  0.05  gm/cm 
for  the  actual  profile.  This  is  probably  an  extreme  case,  since  I0°C  surface  temperature 
errors  are  unusually  large.  Also,  the  effect  is  much  smaller  for  tropical  or  middle 
latitudes,  since  the  total  water  would  be  greater,  and  the  humidity  sensing  channels  are 
then  less  sensitive  to  low-level  temperature  errors. 

3.5  INFORMATION  RETRIEVABLE  FROM  WINDOW  CHANNELS  AT 

795  cm"1  AND  900  cm"1 


The  use  of  window  channel  pairs  to  infer  surface  temperature  and  total 
precipitable  water  implies  that  the  two  measured  radiances  depend,  to  a reasonable 
approximation, on  only  these  two  quantities.  The  approximation  is  expressed  formally  by 
writing  the  radiative  transfer  equation  as 


I.  (u  , T ) — t”.  (u  ) B.  (T  ) + 


N.  (u  ); 
i s’ 


1,2 


(58) 


where  I.  is  the  total  upwelling  radiance  for  channel  i,  B.  (Ts)  is  the  Planck  blackbody  emis- 
sion for  surface  temperature  Ts,  and  t.  (us)  and  Nj  (u$)  are  the  transmittance  and  upward 
radiance,  respectively,  of  the  entire  atmosphere  when  it  contains  total  precipitable  water 
us-  The  last  two  quantities  are  averages  over  a set  of  representative  vertical  water  vapor 
and  temperature  distributions.  The  fact  that  averages  must  be  used  is,  of  course,  the 
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reason  the  equation  is  an  approximation.  The  approximation  becomes  more  accurate  with 

decreasing  u and/or  decreasing  channel  opacity, 
s 

According  to  Equation  58  the  surface  temperature  may  be  expressed  as 


T 

s 


'i  - *1  (US> 
Ti  (us) 


; i = 1,2 


(59) 


where  B.  * (I.)  is  the  inverse  Planck  function.  If  the  radiance  I.  and  empirical  functions 

N.  (u  ),  x'.  (u  ) are  given,  Equation  59  can  be  used  to  obtain  T as  a function  of  u . In  fact, 
ISIS  s s 

two  curves  of  T versus  u can  be  obtained  corresponding  i = I and  i = 2.  The  coordinates  of 

the  intersection  of  the  curves  represent  solutions  for  Tg  and  us*  This  is  essentially  the 

algorithm  used  in  inferring  surface  temperature  and  total  precipitable  water.  When 

Equation  59  is  used  in  this  way,  1.  represents  the  actual  (measured)  radiances,  which  do  not 

necessarily  satisfy  Equation  58  for  the  actual  u and  T ; that  is,  the  solution  values  will 

s s 

generally  be  in  error. 

For  each  of  218  radiosonde  profiles  of  humidity  and  temperature  we  calculated 

the  transmittance  of  the  entire  atmosphere  and  its  radiance  (excluding  the  surface 

contribution)  in  each  of  the  two  window  channels.  Let  these  quantities  be  denoted  t.  , 
k k k 

N,  and  N0  , with  k indicating  the  particular  member  of  the  ensemble  of  218 
£ II  k 

radiosondes.  The  four  quantities  plotted  each  as  a function  u§  , the  total  precipitable 

water  for  ensemble  member  k,  are  shown  as  Figures  70  through  73.  Visually-estimated 

best-fit  curves  were  drawn  through  these  points  to  obtain  the  averages  T.(u  ),  N.(u  ) 

referred  to  above.  These  functions  are  each  represented  in  the  computer  by  a table  of 

65  values  (see  Appendix  A).  Note  that  the  scatter  in  the  plotted  points  becomes  larger  as 

u increases,  corresponding  to  the  deterioration  in  accuracy  of  Equation  58. 
s 

To  test  the  accuracy  of  the  algorithm,  we  used  it  to  retrieve  u and  T for 

5 5 

218x3  cases,  consisting  of  the  218-member  radiosonde  set  combined  with  surface 

k k o k 

temperatures  equal  to  T j and  T | +5  C,  where  T j are  the  radiosonde  air  temperatures 
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just  above  the  surface.  We  realize  that  this  method  of  constructing  surface  temperatures 
is  rather  arbitrary,  but  we  could  not  find  any  data  describing  the  statistics  of  surface 
temperature,  or  possible  correlations  between  effective  surface  temperature  and  low-level 
air  temperature. 

It  was  found  that  the  pair  of  equations  59  do  not  always  have  a unique  solution, 

because  of  the  approximation  inherent  in  using  the  averages  N.  and  t..  That  is,  for  some  of 

the  values  substituted  into  Equation  59  the  two  functions  do  not  intersect,  and  in  some 

cases  there  are  multiple  intersections.  It  was  determined  empirically  that  in  the  case  of 

multiple  intersections,  the  one  corresponding  to  the  largest  u generally  gave  the  best 

solution.  When  there  were  no  intersections  we  determined  u and  T from  the  point  corre- 

s s 

sponding  to  the  minimum  difference  between  the  two  functions. 

Figure  74  shows  the  error  in  retrieved  surface  temperature  as  a function  of  the 

actual  total  precipitable  water  u . Parts  (a),  (b)  and  (c)  of  the  figure  correspond  to  air- 

surface  temperature  discontinuities  of  -5  degrees,  0,  and  +5  degrees  C,  respectively. 

Figure  75  shows  errors  in  retrieved  ug  in  the  same  format.  The  point  symbol  "C"  indicates 

that  Equation  59  had  no  solution  (intersection).  Table  7 gives  the  mean  and  variance  of  the 

errors  as  a function  of  actual  u , and  these  same  statistics  for  the  entire  ensemble;  the 

s 

table  does  not  include  the  cases  for  which  there  was  no  solution. 

Figures  74  and  75  and  Table  7 show  that  the  selected  window  channel  pairs  are 

ineffectual  in  determining  total  precipitable  water:  even  when  the  total  water  is  less  than 

2 o 

I gm/cm  and  the  surface  temperature  discontinuity  is  within  +5°C,  the  mean-plus-one- 

2 

sigma  error  in  inferred  total  water  can  equal  0.5  gm/cm  . However,  the  window  channels 

2 

give  fairly  good  estimates  of  surface  temperature:  for  u <2  gm/cm  the  mean-plus-one- 

o o s 

sigma  error  is  less  than  1C,  and  it  is  less  than  2°C  for  any  of  the  unit  thickness  water  bins 

2 

between  0 and  5 gm/cm  . 

While  these  window  channels  are  not  useful  by  themselves  for  determining  total 

precipitable  water,  they  would  be  excellent  candidates  as  additional  sounding  (inversion) 

channels.  The  weighting  function  for  the  795  cm  * channel  was  plotted  in  the  format  of 

2 

Figure  59,  and  was  found  to  have  a peak  at  a u value  greater  than  4 gm/cm  . 
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Figure  74  a Errors  in  Inferred  Surface  Temperature  Ts  for 
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ERROR  IN  INFERRED  T (C) 


Figure  74  b 
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ERROR  IN  INFERRED  T (C) 


ERROR  IN  INFERRED  o (gm/crn  ) x 10 


Table  7 MEAN  AND  STANDARD  DEVIATIONS  (STD)  OF  ERRORS  IN 
T AND  u INFERRED  FROM  WINDOW  CHANNEL  RADIANCES 

5 3 


u 

s 

0-1 

1-2 

2-3 

3-4 

4-5 

(gm/cm 

AT  = -5°C 

No.  of  Cases* 

35 

38 

33 

41 

20 

171 

Mean  T Error 
$ 

-0.23 

-0.23 

0.05 

0.44 

l.ll 

0.19 

STD  of  Error 

0.14 

0.52 

0.44 

0.62 

0.75 

0.71 

Mean  u$  Error 

0.10 

-0.18 

-0.84 

-1.61 

-0.87 

-0.71 

STD  of  Error 

0.21 

0.74 

0.92 

1.55 

1.50 

1.27 

AT  = 0 

No.  of  Cases* 

35 

32 

27 

32 

19 

148 

Mean  T Error 
s 

-0.22 

-0.24 

-0.22 

0.41 

1.10 

0.09 

STD  of  Error 

0.14 

0.37 

0.72 

1.15 

0.86 

0.86 

Mean  u Error 
s 

0.14 

0.19 

-0.42 

-0.64 

-0.26 

-0.19 

STD  of  Error 

0.22 

1.38 

1.64 

1.33 

0.85 

1.23 

AT  = +5°C 

No.  of  Cases* 

34 

23 

22 

33 

19 

134 

Mean  T Error 

5 

-0.20 

-0.26 

0.61 

0.64 

0.83 

0.27 

STD  of  Error 

0.14 

0.72 

1.12 

1.22 

1.05 

1.04 

Mean  u Error 
$ 

0.21 

-0.01 

-0.01 

-0.20 

-0.65 

-0.09 

STD  of  Error 

0.30 

1.29 

1.46 

0.92 

0.56 

0.99 

Number  of  Cases  for  which  the  Actual  Precipitable  Water  is  in  the  Range  0-1  gm/cm  , 
2 

1-2  gm/cm  , etc.  The  Table  does  not  include  cases  for  which  Eqs.  (59)  had 
no  solution. 
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SECTION  4 

CONCLUSIONS  AND  RECOMMENDATIONS 


4. 1 TEMPERATURE  RETRIEVAL  CAPABILITY 

The  study  compared  the  direct,  nonlinear,  nonstatistical  techniques  of  Smith*, 

2 3 4 

Chow  , Fleming  , a modified  version  of  Twomey's  method  , and  two  variants  of  these 

techniques.  Smith's  iterative  method  was  judged  to  give  the  best  temperature  retrievals, 

for  this  class  of  retrieval  technique. 

3 

Smith's  method  was  then  compared  to  the  Fleming  statistical  direct  method 

and  two  linear  inverse  methods:  the  minimum-information  method  of  Foster^  and 

7 6 8 

Twomey  , and  the  full-statistics  method  of  Foster  and  Strand  and  Westwater  . This  four- 
way comparison  produced  the  following  conclusions: 

1.  The  most  accurate  of  the  four  methods  is  usually  either  the 
full-statistics  method  or  the  Fleming  statistical  method, 
depending  on  the  climatology  represented  by  the  actual 
profile  and  assumptions  made  as  to  the  degree  of  correlation 
between  the  surface  temperature  and  low-level  air  temper- 
ature. From  the  data  available  it  was  not  possible  to 
determine  the  air-surface  temperature  correlation. 

2.  When  sensitivity  to  the  use  of  inappropriate  statistics  and 
ability  to  retrieve  anomalous  profiles  are  included  as  factors 
in  the  evaluation,  the  full-statistics  method  emerges  as  the 
best  retrieval  technique.  We  refer  here  to  statistics  defining 
the  temperature  covariances  between  all  atmospheric  levels, 
not  just  the  air-surface  covariance. 
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3.  Quite  satisfactory  temperature  retrievals  are  obtained  w'th 
either  the  full-statistics  or  Fleming  statistical  methods  using 
only  two  sets  of  statistics,  one  for  the  tropics  and  one  for 
extratropical  latitudes. 

4.  For  the  two  best  methods,  maximum  errors  in  the  temper- 
ature retrievals  below  the  70  mb  level  at  extratropical  lati- 
tudes are  approximately  3 to  4 degrees  C,  when  noise  in  the 
data  is  negligible  (or  effectively  suppressed  by  the  regu- 
larization parameter  values  adopted  in  the  study).  For  the 
tropics  the  maximum  error  below  70  mb  is  4 to  5°C  and 
occurs  near  the  tropopause.  However,  when  the  tropopause  is 
unusually  high  or  low  errors  up  to  8 or  10  degrees  will  occur. 

4.2  HUMIDITY  RETRIEVAL  CAPABILITY 

The  conclusions  relative  to  the  capability  of  DMSP-H  to  retrieve  precipitable 
water  overburden  profiles  are: 

1.  Three  of  the  six  channels  are  essentially  redundant.  The  use 
of  channels  I,  2 and  4,  at  353  cm”*,  347  cm”*  and  408  cm  * 
yields  as  much  information  as  all  six  channels. 

2.  The  channels  are  too  strongly  absorbing  to  infer  low-level 

moisture.  Channel  4,  the  most  transparent,  becomes  opaque 

when  the  precipitable  water  exceeds  approximately 
2 

1 .3  gm/cm  , which  is  a typical  value  for  the  total  water  at 

middle  latitudes.  In  the  tropics  the  total  water  nearly  always 
2 

exceeds  1.3  gm/cm  . 

3.  None  of  the  channels  is  sufficiently  opaque  to  infer  moisture 
above  the  300  mb  level  in  dry  artic  cases. 
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4.  The  full-statistics  method  appears  to  give  better  solutions  for 

the  log  of  the  integrated  water  profile  (In  u versus  p)  than  the 

minimum-information  method.  The  former  method  is  at  least 

able  to  retrieve  In  u values  different  from  the  guess  at 

2 

altitudes  where  the  actual  u is  greater  than  1.3  gm/cm  ; this 
"information”  is  not  really  retrieved  but  comes  from  the 
statistics  represented  by  the  covariance  matrix  and  the  mean 
profile  used  for  the  guess. 

5.  The  In  u profiles  retrieved  by  the  full-statistics  method  retain 
much  of  the  general  shape  characteristics  of  the  guess  In  u 
profile.  Thus,  large  errors  occur  in  the  retrieved  profile  when 
the  actual  profile  differs  significantly  in  shape  from  the  guess 
profile.  However,  the  retrieved  profile  is  unaffected  by 
simple  scaling  of  the  guess  profile.  Errors  in  retrieved  u 
below  the  300  mb  level  were  less  than  30  percent  in  the  best 
of  the  cases  tested,  and  of  the  order  of  100  percent  in  the 
worst  case.  It  is  possible  that  these  errors  can  be  reduced  by 
using  mean  (guess)  profiles  and  covariance  matrices  represen- 
tative of  different  seasons  and  climate  region  as  well  as 
latitude. 

6.  The  relative  humidity  profile  corresponding  to  the  retrieved 
In  u profile  tends  in  many  cases  to  have  a large  oscillatory 
component,  and  bears  little  resemblance  to  the  actual  humid- 
ity profile.  This  is  probably  due  to  the  fact  that  the 
retrievable  information  is  contained  in  only  three  channels. 

7.  The  bands  at  795  cm  * and  900  cm-*  (each  of  width  20  cm  *) 
are  not  sufficiently  transparent  to  be  used  for  accurately 
determining  total  precipitable  water,  although  they  yield 
good  estimates  of  surfoce  temperature  for  middle  and  high 
latitudes  (u<2  gm/cm  ).  These  bands  would  be  good  can- 
didates as  sounding  channels  to  replace  the  redundant  DMSP 
channels. 
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RECOMMENDATIONS 


The  results  of  this  effort  suggest  that  the  following  study  tasks  should  be 
pursued  in  future  research  aimed  at  improving  the  humidity-sensing  performance  of  DMSP: 

1.  Analyze  the  performance  of  DMSP  with  redefined  humidity 
sensing  channels,  consisting  of  the  present  353  cm  * , 

347  cm  ' and  408  cm  ' channels  in  combination  with  bands 
of  width  5 cm  * or  10  cm  * at  795  cm  * and  900  cm  * , and  a 
sixth  band  somewhat  more  opaque  than  the  353  cm  * channel. 

The  weighting  function  -dt/d  In  u*  of  this  sixth  channel 

_4  2 

should  have  a maximum  at  u*  =*  5 x 10  g m/cm  . Consider 
other  proposed  combinations  of  channels  in  the  H^O  rota- 
tional band. 

2.  Analyze  humidity  sensing  performance  when  the  primary 
sensing  channels  are  selected  in  the  6.3  um  water  band  rather 
than  the  rotational  band.  The  weaker  continuum  in  the 
6.3  pm  region  should  result  in  sharper  weighting  functions. 

3.  Investigate  the  use  of  the  two  most  transparent  humidity 
channels  for  determining  surface  temperature.  If  these  bands 
can  determine  surface  temperature  more  accurately  than  the 
six  CC>2  channels,  low-level  errors  in  the  retrieved  temper- 
ature profile  and  corresponding  errors  in  the  retrieved  In  u 
profile  would  be  reduced. 

4.  Study  the  climatological  statistics  of  water  vapor  profiles  to 
determine  whether  distinctly  different  means  and  covariances 
can  be  constructed  according  to  season,  latitude  and  climate 
type.  More  accurate  initial  guesses  and  covariance  matrices 
would  lead  to  improved  retrieval  solutions.  If  sufficient  data 
is  available,  determine  whether  there  is  a significant  correla- 
tion between  surface  temperature  and  low-level  air 
temperature. 
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5.  Study  the  statistics  of  HjO  retrieval  errors,  with  instrument 
noise  effects  included,  and  determine  the  distribution  of 
errors  according  to  climatological  category. 

6.  Study  In  detail  the  relationship  between  temperature  retrie- 
val errors  and  humidity  retrieval  (In  u)  errors. 

7.  Study  the  effects  of  possible  errors  in  the  CC^  and 
transmittance  functions. 

8.  Study  the  effectiveness  of  proposed  techniques  for  correcting 
for  partial  cloud  cover. 
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APPENDIX  A 

DERIVATION  OF  EQUATION  22 


We  seek  the  J x M matrix  C which  best  fits  the  relation 


AB  = CAI  (B  = B ; I = I ) 

^ ^ /v  ^wl  '>*1 

in  the  sense  of  least  squares.  That  is,  C should  minimize  the  error 

J K / M v 2 

««=>■£  - 

j=l  k=l  ' i=l  ' 


(A.I) 


(A  .2) 


in  Eq  A.I  when  it  is  applied  to  the  data  ensembles  AB.^  = B.^  - (Bj>,  Al.^  a 1.^  — The 
columns  of  the  matrices  AB  and  Al  correspond  to  a particular  ensemble  member,  and 
represent  the  departure  from  a mean  Planck  radiance  profile  (at  the  reference  frequency  r) 
and  the  corresponding  departures  from  the  mean  channel  radiances. 

The  solution  for  C is  obtained  from  the  normal  equations  d£/dCnm  = Os 


0 = 


or. 


K , M \ 

KT~  = 2 2 A,mk  lABnk"  2 Cni  Alik) 

nm  k=l  ' U\  ' 

K M , K v 

X]  ABnk  Almk“  X^  Cni(X^  A,ik  Almk) 

k=l  i=l  ' k= I ' 


= AB(AI)T  - C(AIXAI)T  . 


(A. 3) 


► -nr 


Solving  for  C we  obtain 


C = AB(AI)T  [ai(ai)t] 
13 8 
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(A. 4) 


Equations  A.  I and  A .4  represent  a linear  regression  solution  for  AB.  Note  that  the 
radiative  transfer  equation  has  not  been  used  if  both  Al  and  AB  are  measured  data. 

The  vectors  Al  and  AB  are  related  by  the  radiative  transfer  equation,  as  are  the 
data  ensembles  Al  and  AB.  If  we  acknowledge  that  the  elements  Al.^  contain  noise  E^,  the 
relationship  is 


Al  = W(AB)  + E . (A.5) 

Thus, 

C4l<4l>T]im  = L(E  Wij  4Bik  * E,k)(E  % 4B,'k  * Emk) 

= Z^tZfe4^48^)^]  EikEmk  • (A'« 

i j ' k ' k 

from  which 

AI(AI)T  = W[AB(AB)T]WT  + EET  . (A.7) 


Equations  A .6  and  A.7  assume  that  E is  uncorrelated  with  AB,  and  that  the  E.^  are  random 
with  zero  mean.  We  can  show  similarly  that 


AB(A|)T  = [ab(ab)t]wt.  (A  .8) 


According  to  the  second  assumption  stated  above,  EE^  is  (K-I)S^,  where  is 
the  covariance  matrix  of  E.  Also,  AB(AB)T  is  (K-I)S,  where  S is  the  covariance  matrix 
of  B.  Using  these  definitions,  and  substituting  A.7  and  A.8  into  A.4,  we  find  that 


(A.9) 


Equations  A.9  and  A.  I yield  Equation  22  in  Section  2.2.2. 
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APPENDIX  B 


AVERAGE  WEIGHTING  FUNCTIONS  FOR  THE  DMSP  TEMPERATURE 
SENSING  CHANNELS,  AND  EMPIRICAL  FUNCTIONS  N (u^  FOR  THE 
WINDOW  CHANNELS  AT  795  AND  900  cm"1 


Table  B- 1 CHAMMEL  WEIGHTING  FUNCTIONS* 


Channel 


Pressure  I 2 3 4 5 6 


(mb) 

(668) 

(676) 

(695) 

(707) 

(727) 

(747  cm"  ) 

.100 

s. 

0. 

0. 

0. 

0. 

0. 

.154 

.A98IE-02 

.1520E-02 

.6486E-03 

.6295F-04 

* 1 057E-O4 

•1065E-O3 

.225 

.992tE-02 

.1957F-02 

.6519E-03 

•9733E-04 

•3692E-04 

•1339F-03 

.31 7 

.1178E-01 

.2106E-02 

.8407E-03 

* 1 334E-03 

•S376E-04 

•1599E-03 

.434 

• 1»60E-O1 

.1974E-02 

• 111  4E-02 

•2046E-03 

.6673E-04 

•1680E-O3 

*579 

.1382E-01 

•2034F-02 

• 1 508E-02 

•3222E-03 

•9046E-04 

•1694E-03 

.755 

.152IE-01 

.2552E-02 

•2096E-02 

.5141E-03 

•1615E-03 

•1670E-03 

.966 

.1661E-01 

•3289E-02 

.2693E-02 

.742BE-03 

•2531E-03 

•1770E-03 

1.217 

. 1 716E-01 

.4007E-02 

.3117E-02 

.9675F-03 

.3667E-03 

.2028E-03 

1.511 

•1760E-01 

.4701E-02 

.3467E-02 

• 1 1 99E-02 

•5046E-03 

•2316E-03 

1 .852 

.1799E-01 

.5409E-02 

.3677E-02 

• 1 404E-02 

•6751E-03 

.2680E-03 

2.266 

.1916E-01 

.6181E-02 

• 3802E-02 

. 1 581E-02 

.8658E-03 

.3076E-03 

2.696 

.1921E-01 

•7003E-02 

.3903E-02 

. 1 706E-02 

• 1 055E-02 

.3533E-03 

3.206 

.1929E-01 

.7067E-O2 

.4032E-02 

. 1 909E-02 

• 1 247E-02 

.4016E-03 

3.792 

.106tE-Ol 

.8849F-02 

.4256E-02 

.1917E-02 

.1449E-02 

.4585E-03 

6.629 

• 1 980E-0 1 

.9853E-02 

.4524E-02 

•2034E-02 

.1644E-02 

•5112F-03 

5.151 

•IQlOE-Ol 

.1084E-01 

.4814E-02 

•2163E-02 

.1011E-O2 

.5527E-03 

5.956 

.1969E-01 

.1191E-01 

•51 70E-02 

. 2338E-02 

. 1 908E-O2 

.5990E-03 

6.862 

.2022E-01 

.1293E-01 

.5539E-02 

. 2530E-02 

•2154E-02 

•6421E-03 

7.821 

.2083E-O1 

.1409E-01 

.5949E-02 

•2765E-02 

•2325E-02 

.6837E-03 

0.896 

•ei6|E-01 

.1519E-01 

•6357E-02 

•2997E-02 

.2483E-02 

.7155E-03 

18^0^3 

.»190E-01 
•8?28E-0 1 

.1623E-01 

.6744E-02 

•3212E-02 

.2630E-02 

.7396E-03 

11.9*6 

12.753 

.1741E-01 

.7231E-02 

• 3455E-02 

•2807E-02 

.7581E-03 

.*254f-01 

.1859E-01 

.7747E-02 

. 3700E-02 

.2994E-02 

.7727E-03 

16.268 

.2264E-01 

.1980E-01 

.8310E-02 

• 3947E-02 

.3192E-02 

.7043E-O3 

1 V.  9 0 8 

.2268E-01 

•2106F-01 

.8931E-02 

.4212E-02 

•3403E-02 

.7992E-03 

17.678 

•2260E-O1 

.2211E-01 

•9557E-02 

•4451E-02 

.3591E-02 

.8216E-03 

19.586 

.2269E-01 

.2312E-01 

.1020E-01 

•4689E-02 

.3778E-02 

.0486E-O3 

21.633 

•2731E-01 

.2411E-01 

.1090E-01 

.4935E-02 

.3961E-02 

.P816E-03 

23.831 

•2?08E-01 

•2509E-01 

. 1 168E-0 1 

.51 87E-02 

•4143E-0? 

.9223E-03 

26.196 

•2184E-01 

•2601E-01 

.1249E-01 

•5435E-02 

.4316E-02 

.9646E-03 

28.698 

.2153E-01 

.2686E-01 

•1332E-01 

•5665E-02 

.4465E-02 

•1015E-02 

31.390 

• 2i 16E-01 

.2757E-01 

.1416E-01 

•5865E-02 

.4584E-02 

.1069E-02 

36.236 

.2077E-01 

.201OE-O1 

.1492E-01 

•6034E-02 

•4676E-02 

•1122E-02 

37.273 

.2O30E-O1 

.2951E-01 

.1561E-01 

.61 84E-02 

.4755E-02 

•1174E-02 

60.649 

.1Q94E-01 

.2078E-O1 

. 1 624E-0 1 

•6319E-02 

.4R21E-0? 

.1228E-02 

63.918 

.1O53E-01 

.2890F-01 

.1673E-01 

.6411E-02 

.4058E-O2 

•1267E-02 

67.560 

. 1 91SE-01 

.2899E-01 

.1721E-01 

•6505E-02 

.4095E-O2 

•1306E-O2 

51.369 

•] 958E-01 

.2990E-01 

•1766E-01 

•6623E-02 

•4939E-02 

.1353E-02 

55.616 

• 1 794E-01 

.2052E-O1 

. 1801E-01 

•6753E-02 

.4982E-02 

•1400E-02 

59.692 

• 1 f32E-01 

.2804E-01 

.1825E-01 

•6967E-02 

•5015E-02 

•1441E-02 

66.179 

.1669E-01 

.2754E-01 

. 1950E-01 

.6974E-02 

•5044E-02 

•1481E-02 

69.913 

.1600E-01 

.2703E-01 

.1973E-01 

.7080E-02 

.5074E-0? 

.1520E-02 

73.892 

• 153fE“01 

•2624E-01 

. 1900E-01 

•7303E-02 

.5161E-02 

.1595E-02 

79.122 

•146TE-01 

.2530E-01 

. 1 929E-01 

•7581E-02 

.5275E-02 

•1688E-02 

86.612 

.1359E-01 

.2622E-01 

.1962E-01 

.7936E-02 

.5432E-02 

.1803E-O2 

90.369 

.1?7*E-01 

.2315E-01 

. 1 996E-01 

.A291E-02 

.5590E-02 

•1917E-02 

96.399 

.118TE-01 

.2210E-01 

•2029E-01 

.9664E-02 

.5747E-02 

•2031E-O2 

* Entries  are  [(t(p.  . 

L j + l 
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Table  &-J  CHANNEL  WEIGHTING  FUNCTIONS  (Continued) 


Channel 


Pressure 

1 

2 

3 

4 

5 

6 L 

(mb) 

(668) 

(676) 

(695) 

(707) 

(727) 

(747  cm'1) 

10?. 713 

.1 lOlE-Ol 

.2094E-01 

.2065F-01 

.91 35E-02 

•5965E-02 

.2198E-02 

109.315 

•1«13E-01 

. 1 958E  -*01 

.2099E-01 

•9R23E-02 

.6267E-02 

.2442E-02 

115.216 

,9?49E-02 

.1817E-01 

.2133E-01 

• 1 055E-01 

•6585E-02 

.2703E-02 

123.423 

.8374E-02 

.1674E-01 

.2166E-01 

.11 31E-01 

•6913E-0? 

.2974E-02 

130.943 

.7537E-02 

•1537E-01 

.2197E-01 

. 1204E-01 

•7230E-02 

.3242E-02 

138.785 

.6768E-02 

.1409E-01 

•2222E-01 

• 1 277E-01 

• 7535E-02 

•3514E-02 

146.956 

.6O22E-02 

.1282E-01 

.2243E-01 

. 1 349E-0 1 

.7838E-02 

.3791E-02 

155.456 

•5?61E-02 

• 1 1 39F-0 1 

•2249E-01 

•1442E-01 

.8226E-02 

.4173E-02 

184.3?? 

.4526E-0? 

.9897E-02 

.2237E-01 

•1552E-01 

.8678E-02 

.4648E-02 

1 7 3.5-32 

.3810E-0? 

.8429E-02 

•2224E-01 

• 1 661E-01 

.9127E-02 

•5125E-02 

183.106 

.3244E-02 

.7218E-02 

.2196E-01 

. 1 759E-01 

.9652E-02 

.5596E-"2 

193.051 

.2738E-02 

.6099E-02 

.2158E-01 

• 1 855E-01 

.9982E-02 

•6089E-02 

?0  3. 376 

•2313E-02 

.5117E-02 

•211 2E-01 

. 1 949E-0 1 

.1041E-01 

.6603E-O2 

214.089 

.1984E-02 

.4309E-02 

•2050E-01 

•2041E-01 

• 1 085E-01 

.7164E-02 

??5.200 

. 1 702E-02 

.3583E-02 

.1979E-01 

•?129E-01 

. 1 129E-01 

.7734E-02 

236.717 

.1431E-02 

.2883E-02 

.1907E-01 

.2214E-01 

.1173E-D1 

.8298E-02 

?48.648 

.1195E-02 

•2263E-02 

. 1831E-01 

•??89E-01 

.121 4E-i) 1 

.8934E-0? 

261.002 

• 1 022E-02 

.1798E-02 

. 1 739E-01 

.2353E-01 

• 1260E-01 

.9432E-02 

273.789 

.8672E-03 

• 1 373E-02 

. 1 644E-0 1 

•?4l7E-01 

• 1 308E-0 1 

• 1 005E-0 1 

287.018 

.7306E-03 

.9898E-03 

.1546E-01 

•2483E-01 

•1361E-01 

.1072E-01 

300.697 

.6?47E-03 

.6979E-03 

. 1 442E-0 1 

• ?539E-0 1 

.1418E-01 

•1142E-01 

314.935 

.5722E-03 

.5281E-03 

. 1 325E-01 

•?554E-01 

•1474E-01 

.1214E-01 

329. <.*1 

.5346E-03 

.4099E-03 

.1209E-01 

•2567E-01 

.1532E-01 

.1287E-01 

344. 5?6 

.4Q98E-03 

•3033E-03 

. 1 092E-0 1 

.?569E-01 

• 1588E-01 

.1358E-01 

360.097 

.4490E-03 

.2104E-03 

•9786E-02 

.?562E-0I 

•1643E-01 

.1427E-01 

376.165 

.♦415E-03 

•1297E-03 

•8676E-02 

.2545E-01 

• 1 695E-01 

•1495E-01 

392.738 

.4209E-03 

.7409F-04 

.7613E-02 

•?506E-01 

.1742E-01 

.1561E-01 

409. 8?7 

.4O55E-03 

.3749E-04 

•6586E-02 

.?442E-01 

.1787E-01 

•1629E-01 

427.440 

•3974E-03 

.2753E-04 

.5634E-02 

. ?349E-0 1 

.1824E-01 

. 1 694E-0 1 

44$.5R7 

.3894E-03 

.1922E-04 

.4715E-02 

•??5lE-01 

. 1859E-01 

.1757E-01 

464.279 

.3814E-03 

.1176F-04 

•3853E-02 

•2145E-01 

• 1 892E-01 

•1817E-01 

483. 5?3 

.3738E-03 

.5686E-05 

•3062E-02 

•2032E-01 

.1923E-01 

.1873E-01 

503.331 

.3666E-03 

.1487E-05 

.2386E-02 

• 1 91 1E-0 1 

.1953E-01 

. 1 924E-0 1 

523.711 

.3616E-03 

•9076E-06 

.1947E-02 

• 1 785E-01 

.1973E-01 

.1971E-01 

544.675 

.3R69E-03 

.5923E-06 

• 1 545E-02 

. 1 659E-0 1 

. 1 993E-0 1 

.2019E-01 

566.231 

.3626E-03 

.4138F-06 

.1221E-02 

. 1 531E-01 

.2008E-01 

•2065E-01 

588.389 

.3488E-03 

.2744F-06 

.9364E-03 

•1 401E-01 

•2022E-0 1 

.?1 15F-01 

611.160 

.3447E-03 

•1412E-06 

.6679E-03 

. 1 ?70E-0 1 

•2032E-01 

•2162E-01 

634.554 

.3403E-03 

.3823F-07 

.4793E-03 

.1 149E-01 

•2029E-01 

•2196E-01 

658.580 

.3429E-03 

.2584E-07 

•3581E-03 

• 1 034E-01 

.2038E-01 

•2271E-01 

683.250 

.3462E-03 

.1392E-07 

.2505E-03 

.9358E-02 

•2049E-01 

.2348E-01 

708.572 

.3540E-03 

.3968E-08 

. 1 542E-03 

.8212E-02 

.2070E-01 

.2455E-01 

734.558 

.3745E-03 

•1778E-08 

. 1 051E-03 

• 7225E-02 

•2105E-01 

.2633E-01 

761.217 

.4O29E-03 

.9588F-09 

.6641E-04 

•6301E-02 

- .2159E-01 

•2856E-01 

788.561 

•4651E-03 

.4016E-09 

.3693E-04 

.5504E-02 

.2282E-01 

.3234E-01 

816.598 

.5O48E-03 

.9328E-10 

•1917E-04 

•4821E-02 

•2373E-01 

.3570E-01 

845.341 

.5405E-03 

.6257E-10 

.1332E-04 

.4P50E-02 

.2451E-01 

•3907E-01 

874.800 

.6?15E-03 

.3331F-10 

.8115E-05 

•3604E-02 

.2470E-01 

.4175E-01 

904.984 

.6388E-03 

.5292E-1 1 

.3544E-0S 

n .??53E-02 

.2374E-01 

•4158E-01 

935.905 

•6230E-03 

.2577E-11 

.2432E-05 

•2415E-02 

.2190E-01 

.3913E-01 

967.574 

.5858E-03 

•1606E-11 

.1663E-05 

• 1 92BE-02 

•1950E-01 

.3553E-01 

1 000.000 

.5O94E-03 

<9240E-13 

.639BE-06 

.1 339E-02 

.1567E-01 

.2932E-01 

1 000.000 

.21 16E-02 

.1439E-14 

.2237E-06 

.2449E-02 

.4499E-01 

.9664E-01 
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THE  EMPIRICAL  FUNCTIONS  T(u)  AND  N(uc) 


ftllE-05 


